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Abstract 

The  Air  Force  can  save  thousands  of  dollars  by  reducing  the  number  of  blade 
hours  on  the  CH-47  through  finding  an  optimal  mixture  of  CH-47s  and  C-130s  to  conduct 
current  operations  in  Afghanistan  and  Iraq.  Ultimately  these  savings  will  relieve 
maintenance  operations  for  the  CH-47  and  lengthen  the  lifespan  of  the  CH-47  airframe. 
Moreover,  incorporating  C-130s  into  the  operations  will  reduce  cargo  transit  time  from 
supply  depots. 

This  study  looks  at  the  involvement  of  the  C-130  in  CH-47  airlift  operations  to 
reduce  CH-47  usage  and  increase  supply  efficiency.  The  research  focus  is  narrowed  to 
current  airlift  operations  in  Afghanistan  and  Iraq  in  the  CENTCOM  theater  of  operation. 
A  mathematical  representation  of  current  CH-47  operations  augmented  with  C-130s  is 
the  foundation  of  this  research.  Particularly,  these  operations  in  CENTCOM’s  area  of 
operations  are  formulated  as  linear  transportation  problems  using  network  mathematics. 

The  uniqueness  this  research  offers  entails  modified  scenarios  of  the 
transportation  problem  solved  as  an  optimization  model.  AMC  requires  additional 
constraints  to  be  augmented  with  the  basic  transportation  linear  model  that  pushes  this 
application  in  new  areas.  In  addition,  the  uncommon  layout  of  supply  depots  to  the 
specified  receiving  airfields  in  Afghanistan  and  Iraq  provide  an  altogether  new  kind  of 
transportation  problem. 
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OPTIMAL  CH-47  AND  C-130  WORKLOAD  BALANCE 


I.  Introduction 


Background 

The  CH-47  network  inU.S.  Central  Command’s  (CENTCOM’s)  area  of 
operations  is  eomprised  often  airfields.  Two  of  these  airfields  are  eapable  of  handling  C- 
130s.  Additionally,  one  of  these  C-130  eapable  airfields  aets  as  the  eentral  hub  for  the 
area  of  operations  for  Afghanistan  and  Iraq.  Forty  pereent  of  CENTCOM’s  supply 
requirements  for  the  CH-47  move  out  from  the  eentral  hub.  Another  forty  pereent  of 
requirements  are  flown  into  the  eentral  hub  airfield.  The  remaining  twenty  pereent  of 
supplies  and  passengers  are  moved  to  and  from  other  airfields  that  do  not  inelude  the 
eentral  hub.  The  C-130  eapable  airfields  are  utilized  as  standard  CH-47  airfields; 
utilizing  C-130s  to  perform  the  portion  of  CH-47  workload  involving  these  airfields  has 
potential  to  reduee  CH-47  hours. 

Problem  Statement 

The  eurrent  OPS  tempo  and  the  austere  terrain  of  Afghanistan  and  Iraq  plaee  a 
heavy  burden  on  the  CH-47  fleet.  Currently  the  CH-47  sehedule  is  built  independently  of 
the  C-130  sehedule  in  both  Iraq  and  Afghanistan.  Some  of  the  intra-theater  lift 
performed  by  the  CH-47  fleet  is  to  airfields  that  support  C-130  operations.  The  potential 
exists  for  C-130s  to  relieve  some  of  the  CH-47  workload  and  save  blade  hours.  The 
problem  then  beeomes  a  question  of  optimizing  the  blending  of  CH-47  and  C-130  routes. 
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Research  Objectives 

This  research  captures  the  underlying  CH-47  network  structures  in  both  Iraq  and 
Afghanistan.  Additionally,  it  is  desired  to  simulate  flight  times  and  cargo/passenger 
requirements  at  these  locations.  Air  Mobility  Command/ Analysis,  Assessment,  and 
Lessons  Learned  Directorate  (AMC/A9)  can  provide  the  CH-47  and  C-130  flight 
schedules  as  well  as  the  requirements  that  drove  those  schedules.  Once  the  network  is 
defined  along  with  the  requirements  that  drive  the  need  for  CH-47  and  C-130  airlift 
assets,  an  optimization  model  is  developed  to  minimize  costs.  Given  AMC/A9  data  and 
the  optimization  model,  the  following  are  research  objectives: 

a)  Determine  how  many  CH-47  blade  hours  can  be  saved 

b)  Determine  what  additional  C-130  effort  is  required  to  gamer  these  savings 

c)  Build  a  method  to  determine  the  correct  C-130-CH-47  mix  based  on  cargo  and 
passenger  requirements  and  the  number  of  C-130  capable  and  non  C-130 
capable  airfields. 

Research  Focus 

This  research  centers  on  CH-47  transport  requirements  between  the  primary  hub 
and  several  forward  airfields  in  the  CENTCOM  area  of  operations.  Historical  data 
provided  by  Air  Mobility  Command  indicates  the  usual  CH-47  airlift  requirements  since 
the  start  of  the  Afghanistan  and  Iraq  wars  on  terrorism.  This  data  forms  the  foundation  of 
research  analysis  as  the  movement  of  cargo  and  passengers  is  tracked  from  the  central 
hub  to  various  airfields. 
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Methodology. 

The  analysis  is  accomplished  with  the  aid  of  an  AMC  optimization  model  that 
supports  CH-47s  and  C-130s  and  a  developed  Microsoft  Excel®  compatible  optimization 
model.  A  collection  of  flight  schedules  flown  to  Afghanistan  and  Iraq  over  a  specified 
time  period  is  the  data  needed  to  execute  this  research  analysis.  The  essential 
optimization  parameters  include  the  central  and  destination  airfield  locations,  the  number 
of  passenger  requirements,  the  amount  of  cargo  requirements,  and  travel  time  between 
airfields.  AMC/A9  provided  the  bulk  of  the  data  sources. 

Assumptions/Limitations. 

The  developed  optimization  model  is  unique  to  this  transportation  scenario  in 
Afghanistan  and  Iraq.  Moreover,  the  parameters  utilized  are  specific  to  the  CH-47  and 
the  C-130.  It  is  assumed  that  there  are  available  C-130s  to  augment  CH-47  operations  as 
well  as  maintain  a  similar  operation  tempo.  Based  on  the  reduced  number  of  CH-47s 
utilized,  the  number  of  blade  hours  saved  is  estimated  from  the  equivalent  workload 
amount  performed  by  the  C-130s  that  would  otherwise  be  performed  by  CH-47s. 

Implications. 

The  goal  of  this  research  project  is  to  present  Air  Force  and  Army  leadership  with 
an  efficient,  combined  employment  of  joint  airlift  assets.  Results  shape  the  size  of  airlift 
capabilities  needed  to  satisfy  CENTCOM’s  requirements  and  the  allocation  of  those 
capabilities.  CH-47  savings  become  tangible  in  the  form  of  decreased  funding  and 
maintenance  resources,  reduced  airframe  and  blade  wear,  and  retained  CH-47 
availability. 
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Research  Overview 


This  investigative  study  eommenees  with  an  in-depth  review  of  researeh  linear 
formulations  and  heuristie  approaehes  on  spoke-and-hub  networks.  Several  integer 
programming  formulations  are  presented  in  ehronologieal  order  based  on  the  date 
published.  These  formulations  eonsist  of  p-hub  median/eenter  problems,  p-hub  eovering 
problems,  and  eapaeitated  versus  uneapaeitated  hub  loeation  problems.  The  heuristie 
literature  review  inelude  topies  on  the  Genetie  Algorithm,  Tabu  Seareh,  and  Lagrangean 
heuristie.  The  researeh  then  shifts  to  a  review  on  Vehiele-Routing  Problems  (VRP)  and 
provides  several  helpful  formulations  applieable  to  balanee  the  CH-47  workload  with 
support  from  the  C-130.  An  investigation  on  applieation  of  the  above  areas  to  military 
airlift  eoneludes  the  literature  review,  whieh  is  presented  in  Chapter  2. 

In  Chapter  3,  the  insight  garnered  from  the  researeh  studies  and  data  provided  by 
U.S.  Transportation  Command  (USTRANSCOM)  is  used  to  develop  a  methodology  to 
find  a  balanee  in  the  CH-47 ’s  workload.  Capturing  the  CH-47  seenario  and  network  is 
the  first  step  towards  developing  a  mixed  integer  linear  program  in  the  LINGO  ® 
software  environment.  The  program  development  is  based  on  the  eonstraints  established 
by  the  sponsor  and  the  information  gleaned  from  CH-47  sortie  reeords  from  the  Iraq 
theater  from  2009  to  2010. 

Results  and  analysis  presented  in  Chapter  4,  and  eonelusions  presented  in  Chapter 
5  are  based  on  a  demonstrative  seenario  using  aetual  aireraft  parameters,  and  flying  hour 
eosts  from  DoD  reeords.  The  seenario  ineorporates  the  same  number  of  hubs  and 
available  CH-47  aireraft  eurrently  employed  in  the  Iraq  theater.  A  similar  daily  eargo 
requirement  is  generated  and  a  baseline  CH-47  workload  is  eaptured.  Analysis 
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performed  is  handled  by  inerementally  inserting  C-130  aireraft  into  the  seenario  and 
evaluating  the  reduetion  in  needed  CH-47  flying  hours,  and  impaet  on  daily  eost. 
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II.  Literature  Review 


Introduction 

The  geographical  scenario  presented  by  the  research  problem  at  hand  can 
unquestionably  be  represented  as  a  “hub-and-spoke”  network  design.  This  design  is  very 
common  in  multiple  industries  such  as  network  communications,  package  delivery,  and 
passenger  airlines.  Despite  the  common  network  design,  these  hub-and-spoke  problems 
belong  to  the  NP-hard  class  of  problems  (Abdinnour-Helm  and  Venkataramanan,  1998). 
In  order  to  tackle  these  difficult  problems,  this  literature  review  encompasses  previous 
research  on  solving  hub-location  problems  looking  at  heuristic  and  integer  programming 
approaches.  Additionally,  a  close  look  at  vehicle-routing  problems  is  included  in  the 
following  review  as  they  are  closely  related,  if  not  all  encompassing,  to  hub-location 
problems.  Both  problems  types  can  be,  and  indeed  have  been,  used  to  determine  how  to 
optimize  vehicle  coverage,  but  the  key  difference  is  hub-location  problems  are  aimed  to 
optimize  flow  throughout  a  network  whereas  vehicle-routing  problems  tend  to  optimize 
coverage.  This  research  optimizes  both  passenger  and  cargo  flow  through  the  network  as 
well  as  ensures  that  there  are  enough  air  assets  to  cover  the  locations  that  have 
requirements.  The  last  portion  of  this  literature  review  covers  military  application  and 
research  of  solutions  to  hub-and-spoke  network  designs. 

Hub-and-spoke  Research 

This  part  of  the  review  looks  at  several  heuristic  and  mathematical  programming 
methods  for  solving  hub-location  problems.  Generally,  the  solution  techniques 
developed  can  be  applied  to  scenarios  containing  multiple  hubs.  Some  research 
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discussed  places  emphasis  on  uneapacitated  versus  eapaeitated  hubs  sueh  that  the  flow 
may  be  limited  in  the  latter  ease.  It  is  important  to  note  that  the  bulk  of  the  researeh 
reviewed  treats  the  loeation  and  number  of  hubs  as  variables  whieh  signifieantly 
inereases  the  diffleulty  in  the  problem  to  be  solved.  However,  the  Afghanistan  and  Iraq 
hubs  unique  to  this  researeh  problem  are  known  and  fixed.  This  leads  to  the  researeh 
questions  that  need  to  be  answered; 

•  Whieh  spokes  are  linked  to  eaeh  hub? 

•  How  mueh  of  the  requirements  ean  be  flown  out  of  and  between  the  hubs 
(eapacities)? 

The  approaehes  previously  researehed  and  developed  ean  be  utilized  to  help  answer  these 
questions. 

Integer  Programming  Formulations. 

This  seetion  of  the  literature  review  provides  the  fundamental  baekground  related 
to  the  intended  methodology  to  be  implemented  in  formulating  the  CH-47  hub-and-spoke 
networks  in  Afghanistan  and  Iraq.  The  takeaways  from  the  ehosen  literature  below  are; 
how  to  formulate  a  generie  hub-and-spoke  network  as  an  integer,  preferably,  linear 
program;  how  to  aeeount  for  flow  requirements;  and  the  assignment  of  resourees  needed 
to  deliver  the  requirements  from  hub  to  destination.  This  last  aspeet  is  partieularly 
signifieant  to  the  researeh  as  the  goal  is  to  reduce  CH-47  blade  hours  by  ineorporating 
other  resourees  (C-130s)  to  deliver  requirements.  Onee  the  network  is  eaptured  and 
formulated  as  a  linear  program,  the  only  aspeets  which  change  are  the  daily  flow 
requirements  and  resourees  required. 
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The  first  model  diseussed  is  O’Kelly’s  quadratie  integer  program  through  the 
perspeetive  of  Aykin.  While  this  model  is  formulated  as  a  quadratie  program,  it  sheds 
light  on  eoneepts  that  ean  be  readily  applied  to  the  linear  hub  loeation  program.  The 
general  problem  is  to  find  the  loeations  of  the  hub  faeilities  and  the  node  assignments  that 
minimize  the  total  transportation  eost  (Aykin,  1990).  Again  the  researeh  in  this  thesis  is 
not  interested  in  determining  hub  loeations  as  they  are  fixed  and  known.  What  is  of 
interest  is  the  aireraft  assignment  methodology  for  the  hubs.  In  relation  to  the  researeh,  it 
is  desired  to  find  an  assignment  formulation  for  multiple  airlift  resourees.  This 
formulation  resembles  the  transportation  problem  formulation. 

Variables  of  interest  in  O’Kelly’s  quadratie  integer  program  are  the  hub  loeations 
and  whieh  nodes  are  assigned  to  them.  Given  n  interaeting  nodes,  flows  between  pairs  of 
nodes  denoted  by  Wy,  the  transportation  eost  Cij  of  a  unit  of  flow  between  nodes  i  and  j, 
and  the  number  of  hub  faeilities  p  to  be  loeated  (Aykin,  1990),  the  following  model  is 
developed  to  solve  for  node  assignment  variable  Ai: 


Minimize  Z  =  Sill  (j^ik  Qwi)] 

i  j  k  m 

subject  to  I  Aifc  <  (n  -  p  +  l)Afcfc  for  all  k, 

i 

^  Xii^  =  1  for  all  i, 

k 

^  ^kk  —  P> 

k 

0  <  Xij^  <  1  and  integer  for  all  i  and  k 


Figure  1  -  O’Kelly’s  formulated  model  (Aykin,  1990) 


In  the  above  formulation,  k  represents  the  loeation  of  hub,  but  a  was  not  speoifieally 
defined. 

Campbell  (1994)  presents  four  additional  integer  programming  formulations 
fitting  this  problem.  Two  of  these  formulations,  the  hub  median  and  uneapaeitated  hub 
loeation  problems,  have  been  eited  and  applied  by  researehers.  The  last  two  formulations 
introdueed  by  Campbell  model  the  diserete  hub  center  and  hub  covering  problems.  In 
regards  to  discrete  hub  location  problems,  the  following  five  items  are  considered  to  be 
given,  or  known; 

1)  n  demand  locations 

2)  r  potential  hub  locations 

3)  The  flow  for  the  demand  location  pairs 

4)  The  per  unit  cost  between  all  location  pairs,  and 

5)  The  hub-to-hub  discount  factor  a  (Campbell,  1994). 

The  first  hub-and-spoke  model  discussed  is  the  p-\mb  median  problem.  Its 
objective  is  to  minimize  total  cost  of  flow  units  and  is  considered  to  be  particularly 
important  to  the  airline  industry.  Before  the  model  is  shown  below,  the  following 
variables  and  parameters  have  been  defined  for  the  set  of  hub  location  problems:  Xykm  is 
the  fraction  of  flow  from  location  i  to  location  j  that  is  routed  via  hubs  at  locations  k  and 
m  in  that  order,  Yk  is  a  (1,0)  variable  if  location  k  is  a  hub,  Zjk  is  a  (1,0)  variable  if  location 
i  is  allocated  to  the  hub  at  location  k,  Wy  is  the  flow  from  location  i  to  location  j,  and  Cy  is 
the  standard  cost  per  unit  from  location  i  to  j  (Campbell,  1994).  Given  these  definitions, 
the  below  p-hwb  median  problem  is  formulated  as  a  linear  program  below; 
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Minimize 


1111 

i  j  k  m 

subject  toJ^Y,=  p. 


j^i  jkm  Q  jkm 


k 

0  <  Vfc  <  1  and  integer  for  all  k, 
0  <  Xiji^^  <  1  for  all  i,j,  k,  m, 


^ijkm  =  1  for  all  i,j, 


^ijkm  ^  Yk  for  all  i,j,  k,  m, 

^ijkm  ^  Y^  for  all  i,j,  k,  m 

Figure  2  -  Campbell’s  (1994)  formulated /?-hub  median  problem 
The  next  formulation  is  that  of  the  uneapaeitated  hub  loeation  problem.  The 
underlying  differenee  here  is  that  the  number  of  hubs  is  an  unknown.  Consequently, 
sinee  the  number  of  hubs  is  unknown,  there  is  an  assoeiated  fixed  eost  for  generating 
eaeh  hub.  As  sueh,  Campbell  defines  the  parameter  Fk  as  the  fixed  eost  of  establishing  a 
faeility  at  loeation  k  (Campbell,  2004).  This  model  is  formulated  as  shown  below: 

+  Y,fA 

i  j  k  m  k 

subject  to  0  <  Vfc  <  1  and  integer  for  all  k, 

0  <  Xiji^^  <  1  for  all  i,J,  k,  m, 

^ijkm  =  1  for  all  i,j, 

k  m 

Xijkm  ^  Yj^for  all  i,j,k,m, 

Xijkm  ^  Y^  for  all  i,j,  k,  m 


Minimize  1111  j^i  jkm  Q  jkm 
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Figure  3  -  Campbell’s  (1994)  formulated  uneapaeitated  hub  loeation  problem 
For  the  hub  eenter  problem,  the  goal  is  to  minimize  the  maximum  eost  for  any 
origin-destination  pair  of  nodes.  The  hub  eenter  itself  ean  be  a  single  hub  or  eonsist  of  a 
eolleetion  of  hubs.  This  type  of  hub  eenter  is  important  for  a  hub  involving  perishable  or 
time  sensitive  items,  in  whieh  eost  refers  to  time,  a  is  a  time  diseount  faetor  due  to  higher 
speed  on  the  inter-hub  links,  and  the  maximum  time  from  an  origin-to-destination  is  of 
interest  (Campbell,  1994).  His  basie  /?-hub  eenter  problem  formulation  is  presented 
below; 


.  .  Maximum  r 

MiniTTllze  X^ijkm^ijkmj 


subject  p, 


0  <Yi^  <  1  and  integer  for  all  k, 
=  ^for  all  i,j, 


k  m 


^ijkm  ^  Ykfor  all  i,j,k,m, 

Xijkm  ^  Ym  for  all  i,j,k,m, 

0  <  Xiji^^  <  1  for  all  i,j,  k,m 
Figure  4  -  Campbell’s  (1994)  formulated p-huh  eenter  problem 
Additionally,  several  variants  of  this  formulation  were  developed  by  Campbell.  These 
formulations  ean;  minimize  the  maximum  eost  for  movement  on  any  single  link,  aeeount 
for  flow  thresholds  for  spoke  nodes,  and  take  the  form  of  a  quadratie  program  (Campbell, 
1994). 
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The  last  hub  location  problem  type  developed  is  the  hub  covering  problems, 
which  Campbell  describes  as  being  analogous  to  facility  covering  problems.  It  is  desired 
to  determine  locations  for  the  hubs  such  that  they  can  satisfy  all  destination  node 
demands  while  minimizing  the  cost  of  the  hubs.  In  this  formulation,  the  objective 
function  is  merely  the  hub  cost  piece  from  the  uncapacitated  hub  location  formulation. 


Minimize 


FA 


k 

subject  to  0  <  Tfc  <  1  and  integer  for  all  k 
^ijkm  ^  Yk  for  all  i,j,  k,  m, 

Akm  ^  Y^  for  all  i,j,k,m, 


V  y 


>  1  for  all  i,J 


k  m 

Figure  5  -  Campbell’s  (1994)  formulated  hub  covering  problem 
Aykin  goes  on  to  confirm  the  applicability  of  the  hub-and-spoke  network  design 
to  scenarios  such  as  those  posed  by  the  CH-47  blade  wear  dilemma  by  indicating  that  one 
of  the  successful  applications  of  the  hubbing  concept  is  found  in  air  passenger/cargo 
transportation  (Aykin,  1995).  He  examines  network  designs  in  which  it  is  possible  to 
flow  units  in  a  network  from  demand  point  to  another  demand  point  (i.e.  nonstop  travel) 
or  where  flow  is  channeled  through  hubs  before  reaching  the  final  demand  point 
(transient  travel).  The  goal  of  these  problems  is  to  minimize  total  transportation  costs 
through  the  network  by  determining  an  optimal  mix  of  service  types  and  hub  locations 
(Aykin,  1995).  The  general  hub  location  and  routing  problem  formulated  as  a 
mathematical  model  accounts  for  an  unknown  number  of  hubs  and  whether  direct  or 
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indirect  flow  is  utilized  in  the  network.  The  following  variables  and  parameters  for  such 
a  mode  are  defined  as; 

•  Xij  is  a  binary  variable  in  which  one  is  assigned  if  flows  from  i  to  j  are  shipped 
directly 

•  Xiktj  is  a  binary  variable  in  which  one  is  assigned  if  flows  from  i  to  j  are 
shipped  with  the  routing  i  to  hub  k  to  hub  t  to  j 

•  Wij  is  the  amount  of  flows  from  demand  point  i  to  demand  point  j 

•  Oij  is  a  binary  variable  in  which  one  is  assigned  if  the  pair  {i,j)  is  an  element 
within  the  set  of  routes  where  nonstop  service  is  permitted 

•  d{?,  Q)  is  the  distance  between  nodes  P  and  Q 

•  and  Cij  is  the  associated  cost  parameter  (Aykin,  1995) 

Given  these  definitions,  Aykin’s  general  model  formulation  of  the  problem  is  shown 
below; 


Minimize  11  0,jW,jC,jd{P„Pj)X,j 
i  j 

^1111  Wij  (^a^Cikd{Pi,  Qk)  +  aCktd{Qk,  Qt)  +  a2Ctjd{Qt,Pj)^X, 
j 

subject  to  Oij  Xij  +  11  ^iktj  =  1  for  all  i,j, 


i  k  t  j 


k  t 

Q]^  G  for  all  k,Xij,Xip,t-j  G  [0,1]  and  integer  for  all  i,  k,  t,J 
Figure  6  -  Aykin’s  (1995)  general  hub  location  and  routing  problem 
The  author  goes  on  to  indicate  that  for  known  hub  locations,  the  above  formulation  can 
be  reduced  to  a  decision  model  that  only  accounts  for  service  types; 
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Minimize 


i  j 


1111 


j  ^iktj^ikt  j 


subject  to  OijXij  + 


^iktj  =  1  for  all  i,J, 


k  t 

Xij  and  Xii^f-j  G  [0,1]  and  integer  for  all  i,  k,  t,J 
Figure  7  -  Aykin’s  (1995)  formulation  for  known  hub  location  and  routing  problem 
Five  methods  are  proposed  in  order  to  find  a  starting  solution  to  these  problems. 
The  first  is  to  randomly  choose  initial  hub  locations.  Another  method  is  to  randomly 
select  a  number  of  demand  nodes  to  act  as  the  beginning  hub  locations.  As  a  third 
method,  Aykin  describes  a  ‘Drop  Solution’  in  which  every  demand  node  is  assumed  to  be 
a  hub.  Each  hub  is  then  evaluated  to  determine  the  increase  in  the  objective  function 
value  if  it  is  closed;  this  is  done  by  removing  the  hub  locations  one  at  a  time  (Aykin, 
1995).  Additionally,  he  proposes  a  drop  and  interchange  method  which  incorporates 
obtaining  an  initial  solution  from  the  drop  technique,  but  then  proceeds  to  interchange  the 
initially  selected  hubs  with  non-hub  demand  nodes  to  see  if  better  objective  function 
values  can  be  obtained.  The  last  suggested  method  is  enumeration,  but  it  is  highly 
discouraged  for  use  in  expansive  networks. 

Ernst  and  Krishnamoorthy  provide  a  similar,  but  altered  formulation  of  the  p-hwb 
median  problem  that  was  presented  earlier  in  Campbell’s  research  and  model 
development.  The  primary  alteration  lies  in  the  objective  function  formulation. 
Campbell’s  is  comprised  of  a  single  piece  which  accounts  for  flow  and  cost  between 
nodes.  Ernst’s  and  Krishnammorthy’s  formulation  consists  of  two  distinct  parts:  the  first 
piece  captures  the  flow  and  costs  from  hub  to  demand  node,  and  vice-versa,  and  the 
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second  piece  accounts  for  transfer  costs  between  hubs,  as  applicable.  For  an 
uncapacitated  single  allocation  problem,  the  authors  present  the  following  mixed  integer 
linear  program; 


iEN  kEN 


iEN  kEN  lEN 


Min  11  dikZikOcOi  +  SDi)  +  111 

i 

s.t.^Zfcfc  =  p, 


kEN 


=  l,Vi  G  N, 

kEN 

Zik  <  Zkk.  Vi,  kEN, 

^  ^  Yi  =  OiZik  -  ^  WijZij  Vi,  kEN, 


lEN  lEN 


jEN 

Zik  ^  {0.1}  vi,/c  G  N, 


Yii  >0\/i,k,lEN 

Figure  8  -  Ernst  and  Krishnmoorthy  (1998)  formulated  p-hwb  median  problem 
They  also  present  an  uncapacitated  multiple  allocation  p-hwb  median  problem 
formulation  where  the  model  accounts  for  varying  flow  depending  on  the  direction 
between  two  nodes  (to  include  hubs  as  well).  This  particular  formulation  could  be  very 
useful  in  developing  the  unique  flow  constraints  relative  to  the  CH-47  problem  as  cargo 
and  PAX  requirements  will  vary  depending  on  the  air  assets  being  employed.  Ernst’s  and 
Krishnamoorthy’s  formulation  is; 


Min 


11  Xd-ikZik  T  adkiY^i  +  zz  SdijXlj 

iEN  kEN 


kEN  IEN 


lEN  jEN 
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kEN 


Hk  =  V> 


kEN 


=  Oi,\/i  G  N, 


Y^Xii  =  Wij,Vi,jEN, 

kEN 


lEN  jEN  lEN 


~  ^ik 


=  0,  Vi,  kEN, 


Zik  <  OiH^,  Vi,  kEN, 

Xlj  <  WijHi,yi,i,lEN, 

Zlj.Yki’Zik  ^  0,Vi,;,/c,/  G  N, 

Hk  G  {0,1}  VkEN 

Figure  9  -  Ernst  and  Krishnmoorthy  (1998)  formulated  uneapaeitated  multiple  alloeation 

p-\mb  median  problem 

Reeent  researeh  papers  generally  utilized  the  above  hub-and-spoke  formulations 
as  a  foundation  for  further  model  variants.  For  example,  de  Camargo  et  al.  introdueed  a 
model  formulation  for  multiple  alloeation  hub  loeation  problems  that  aeeounts  for  flow 
eongestion.  This  model  utilizes  a  foundation  formulation  that  largely  resembles  the 
researeh  of  Campbell  and  Ernst  and  Krishnamoorthy.  In  de  Camargo ’s  researeh,  the 
team  defines  a  funetion  that  assesses  hub  loading  and  distributes  the  load  among  other 
hubs  mitigating  flow  eongestion  (de  Camargo,  Miranda,  Ferreira,  Euna,  2009).  Hyun  and 
O’ Kelly  eonstrueted  the  reliable  /?-hub  loeation  problem  that  foeuses  on  loeating  /?-hubs 
on  a  network  to  improve  network  reliability  to  deliver  interaeting  flows  among  its  set  of 
origin-destination  nodes  (Hyun,  O’Kelly,  2009).  The  development  of  their  model  relied 
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on  similar  concepts  for  the  single  and  multiple  alloeation  /?-hub  loeation  problems.  A 
signifieant  eontribution  to  the  set  of  p-hwb  loeation  problems  from  their  researeh  is  the 
inelusion  of  reliability  of  routes  in  the  network.  These  papers  do  not  neeessarily  provide 
relevant  model  formulations  for  the  CH-47  researeh  seenario  of  this  thesis,  but  the 
formulation  insight  and  eoneept  eontributions  are  invaluable  to  methodology 
development. 

Heuristic  Approaches. 

Abdinnour-Helm’s  and  Venkataramanan’s  researeh  into  solution  teehniques  for 
hub-and-spoke  similar  problems  diseovered  that  Morton  O’Kelly  initially  developed  a 
quadratie  integer  program  to  solve  the  p-Yhxh  Median  Problem  (Abdinnour-Helm  and 
Venkataramanan,  1998).  The  heuristies  developed  by  O’Kelly  to  solve  his  integer 
program  involve  eomplete  enumeration  of  all  possible  eonfigurations.  In  his  first 
heuristie,  the  enumeration  proeess  will  yield  an  upper  bound  on  the  objeetive  funetion 
under  the  assumption  that  eaeh  spoke  node  is  attaehed  to  the  nearest  hub  (O’Kelly,  1987). 
The  seeond  heuristie  examines  all  nodes  with  respeet  to  its  first  or  seeond  nearest  hub. 
These  heuristies  offer  important  insight  into  the  problem  at  hand,  but  an  additional 
assumption  must  be  faetored  in  to  fix  the  hubs  in-plaee  and  address  the  alloeation  of 
nodes  to  those  hubs. 

Alternative  approaehes  proposed  by  Klineewiez  involve  exehange  heuristies 
whieh  seek  to  eombine  multiple  metries  to  weed  out  hubs  that  offer  less  performanee 
eompared  to  those  already  found  (Klineewiez,  1991).  The  first  approaeh  is  a  single- 
exehange  heuristie  whieh  orders  a  set  of  spoke  nodes  based  upon  distanee  and  flow  traffie 
and  then  evaluates  those  nodes  with  respeet  to  the  its  nearest  hub(s).  The  new  hub,  onee 
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evaluated  to  a  partieular  new  node  that  offers  performanee  gains,  are  then  swapped  for 
the  old  hub.  The  seeond  approaeh  is  very  similar  to  the  first  exeept  it  examines  a  pair  of 
spoke  nodes  at  a  time  to  new  hubs;  henee,  it  is  ealled  a  double-exehange  heuristie. 
Klineewiez’s  final  proposed  heuristie  is  based  on  elustering.  This  teehnique  arbitrarily 
takes  a  group  (eluster)  of  nodes  and  then  finds  the  best  suited  hub  based  upon  the 
eolleetive  metries  of  the  eluster  (Klineewiez,  1991).  Again,  these  heuristies  are 
developed  under  the  seenario  that  the  hub  loeations  are  unknown.  But  elustering 
teehniques  may  prove  useful  in  solution  development  to  USTRANSCOM’s  goal  of 
redueing  wear  on  CH-47  blades. 

Another  author  goes  on  to  distinguish  between  multiple  alloeation  p-hwb  median 
problems  and  single  alloeation  ones.  Campbell  designates  those  that  handle  multiple 
alloeations  as  HMP  and  those  of  single  alloeation  as  HMP-S  (Campbell,  1995).  Unless 
this  thesis  speeifies  otherwise,  all  p-hwb  median  problems  and  their  solutions  are 
eonsidered  to  be  of  the  multiple  alloeation  type.  In  developing  solution  heuristies  of 
HMP-S,  Campbell  utilizes  the  solution  of  HMP  as  a  beginning  point  in  obtaining  a 
solution  for  HMP-S  (Campbell,  1995).  A  greedy-interehange  heuristie  is  used  to  get  a 
feasible  solution  to  HMP  through  enumeration  and  then  interehanges  hubs  that  result  in  a 
lower  transportation  eost.  Onee  a  solution  to  HMP  is  found,  all  multiple  alloeations  for 
demand  points  must  be  replaeed  by  single  alloeations  (Campbell,  1995).  At  this  point, 
two  heuristies  were  developed  to  obtain  a  solution  to  HMP-S.  The  first  evaluates 
maximum  flow  by  linking  nodes  to  the  loeation  of  the  hub  in  whieh  maximum  flow  is 
generated.  The  seeond  heuristie  minimizes  transportation  eost  by  evaluating  all 
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combinations  of  the  nodes  available.  It  is  important  to  note  that  there  are  no  bounding 
restrietions  on  the  number  of  demand  points. 

Aside  from  p-hwb  median  loeation  problems,  there  is  another  set  of  similar  hub- 
and-spoke  type  problems  ealled  uneapaeitiated  hub  loeation  problems,  or  UHP.  The 
determining  distinetion  of  these  problems  is  that  the  number  of  hubs  is  itself  a  variable. 
Although  the  seenario  of  this  researeh  deals  with  a  known  number  of  hubs  with  known 
loeations,  the  solution  approaehes  to  UHP  are  still  relevant. 

In  order  to  solve  UHP,  Abdinnour  generated  a  new  heuristie  that  blends  a  Genetie 
Algorithm  (GA)  and  Tabu  Seareh  (TS).  Before  any  further  diseussion  is  made  on 
Abdinnour’s  new  heuristie,  whieh  is  ealled  GATS  for  “Genetie  Algorithm  and  Tabu 
Seareh”  (Abdinnour,  1998),  the  following  paragraphs  will  summarize  GA  and  TS. 

GA  is  derived  from  the  idea  of  natural  seleetion  and  geneties  to  produee  a  more 
“evolved”  set  of  solutions  given  an  initial  starting  solution(s).  Eaeh  initial  solution  is 
evaluated  to  determine  how  well  it  performs.  From  here  a  new  set  of  solutions  is 
developed  from  eombining  individual  solutions  from  the  initial  set.  This  eyele  eontinues 
until  a  desired,  or  “no  better,”  solution  is  obtained. 

In  TS,  an  iterative  seareh  proeedure  is  employed  sueh  that  eaeh  iteration  moves 
from  one  feasible  solution  to  another.  Onee  this  move  is  eompleted,  the  algorithm  is 
prohibited  to  ever  go  baek  to  a  previous  feasible  solution  until  a  speeified  number  of 
iterations.  This  is  to  prevent  eyeling.  However,  eaeh  move  may  not  neeessarily  be  an 
improvement  over  the  previous  solution.  The  idea  is  to  avoid  getting  trapped  in  loeal 
optimal  solutions,  i.e.  you  are  foreed  to  test  around  the  feasible  region  for  potentially 
better  solutions. 
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Abdinnour  outlines  the  new  heuristie  in  five  steps.  The  starting  step  involves 
inputting  the  distanee  and  flow  data  for  the  interaeting  nodes  in  the  network,  and 
initializing  the  best  solution  of  GATS  to  a  very  large  number  (Abdinnour,  1998).  The 
next  step  invokes  the  GA  heuristie.  Onee  the  best  solution  is  found,  it  is  saved  and 
passed  on  to  the  third  step,  whieh  is  the  TS  phase.  TS  jumps  around  the  feasible  region 
for  a  predetermined  number  of  iterations,  at  which  point  the  best  solution  is  saved  and 
distinguished  from  the  GA  solution.  The  fourth  step  simply  compares  solutions  between 
the  two  sub-heuristic  phases.  This  comparison  results  in  the  overall  best  GATS  solution 
being  updated  and  then  steps  two  through  four  are  repeated  for  a  predefined  number  of 
iterations.  The  final  step  of  the  GATS  heuristic  is  the  reporting  phase. 

The  last  heuristic  discussed  is  a  Lagrangean  heuristic  developed  by  Elhedhli  and 
Wu.  Their  approach  is  intended  to  solve  hub-and-spoke  networks  formulated  as  a 
nonlinear  mixed-integer  program  that  explicitly  minimizes  congestion,  capacity 
acquisition,  and  transportation  costs  (Elhedhli  and  Wu,  2010).  The  formulated  problem 
is  broken  into  a  more  solvable  subproblem  and  an  NP-Hard  nonlinear  subproblem.  The 
nonlinear  subproblem  is  first  linearized  using  piecewise  functions  and  then  solved  to 
optimality  using  a  cutting  plane  method  (Elhedhli  and  Wu,  2010).  From  the  subproblem 
solutions,  a  heuristic  solution  can  be  found  to  the  original  problem  formulation.  This 
research  by  Elhedhli  and  Wu  is  intended  for  large-scale  networks  and  takes  a  hard  look  at 
mitigating  congestion.  The  scenario  under  research  in  this  thesis  is  not  so  much 
concerned  with  congestion  as  the  underlying  network  is  substantially  large  utilizing  few 
vehicles,  but  it  does  offer  useful  insight  into  hub-and-spoke  problems  modeled  as 
nonlinear  programs. 
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Vehicle-Routing  Problems 

In  the  most  general  sense,  the  vehiele  routing  problem  ean  be  defined  as  a  set  of 
eustomers  eaeh  with  a  known  loeation  and  a  known  requirement  for  some  eommodity, 
whieh  is  to  be  supplied  from  a  set  of  depots  by  a  set  of  delivery  vehieles  of  known 
eapaeity  (Kulkarni  and  Bhave,  1985).  This  definition  exaetly  eaptures  the  overall  goal  of 
this  thesis.  Under  the  hub-and-spoke  design,  our  spoke  (destination)  nodes  are  the 
eustomers  that  have  PAX  and  eargo  requirements  from  a  set  of  hubs.  Currently,  the  CH- 
47  is  the  only  delivery  vehiele  used  in  this  proeess.  The  desire  is  to  add  the  C-130,  and 
potentially  other  eargo  aireraft,  to  the  vehiele  set. 

Kulkarni  and  Bhave  show  that,  under  eertain  assumptions,  the  vehiele-routing 
problem  (VRP)  ean  be  redueed  to  the  travelling  salesman  problem  (TSP).  Using  the  TSP 
as  a  eore  model,  they  show  how  two  formulations  for  the  VRP  were  developed;  single 
and  multiple  depot  VRPs.  In  the  following  formulation,  V represents  the  number  of 
vehieles,  Pk  is  the  eapaeity  of  vehiele  k,  77  is  the  maximum  eost  allowed  for  a  route  of 
vehiele  k,  Qi  is  the  demand  at  node  /  {Qf^  =  0),  and  Xijk  is  a  binary  variable  where  it  is  a 
one  if  the  pair  {i,j)  is  in  the  route  of  vehiele  k  (Kulkarni  and  Bhave,  1985).  Sinee  the 
formulation  for  multiple  depot  VRPs  ean  be  applied  to  seenarios  involving  a  single  depot, 
only  the  multiple  depot  formulation  is  presented  here; 
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_/=W+l  i  =  l 

^ijk  —  0  or  1  for  all  i,j,  k, 

yi  —  Vj  +  (M  +  N)Xiji^  <  N  +  M  —  1  for  1  <  i  ^  j  <  N  and  1  <  k  <V 
Figure  10  -  Kulkarni  and  Bhave  (1985)  formulated  multiple  depot  VRP 
This  formulation  assumes  that  a  single  vehiele  can  satisfy  a  demand  node.  Additionally, 
the  constraints  only  permit  one  vehicle  to  be  assigned  to  a  demand  node.  Kulkarni  and 
Bhave  present  a  new  formulation  of  the  VRP  based  upon  the  above  model  initially 
formulated  by  B.  L.  Golden.  In  their  model,  the  authors  consider  the  case  where  vehicle 
capacity  and  maximum  route  cost  are  identical  per  vehicle.  Kulkarni’s  and  Bhave’s  new 
VRP  formulation  is  given  below; 
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Minimize  Z 


N+M 


N+MN+M 

11 

t=i  i=i 


subject  to  ^  Xij  =  1  for  j  =  1,2, 


1=1 


N+M 


^  Xij  =  1  for  i  =  1,2, 

W+M  W+M 

I  1^0=1'. 

t=W+l  j  =  l 
N+M  N+M 

1  1+J  =  ''' 


j=N+l  i  =  l 


Xij  =  0  or  1  for  all  i,j, 
yi  ~  yj  +  L,^ij  ^  b  —  1  for  1  <  i  ^  j  <  N, 


Ui  —  Uj  +  PXij  <  P  —  Qi  for  1  <  i  ^  J  <  N, 

Vi  —  Vj  +  T  Xij  <T  —  Cij  for  1  <  i  ^  j  <  N, 

Figure  1 1  -  Kulkarni  and  Bhave  (1985)  multiple  depot  VRP  model  for  similar  vehieles 
Researeh  into  the  Defense  Courier  Serviee  distribution  network  by  Baker  (1991) 
explored  the  effeets  of  the  network  by  redueing  the  number  of  depots,  shifting  the 
loeations,  and  ehanging  the  flight  routes  by  formulating  the  network  as  a  VRP  with  an 
underlying  hub-and-spoke  strueture.  The  model  developed  is  largely  based  upon 
Kulkarni ’s  formulation  above  along  with  a  formulation  by  LaPorte  (not  ineluded  in  this 
review).  Baker’s  researeh  led  to  a  model  eapable  of  eapturing  extensive  networks  that 
permitted  multiple  depots  and  tours,  but  due  to  the  expansiveness  of  the  problem,  it  eould 


23 


not  be  solved  exaetly  (Baker,  1991).  In  order  to  find  a  solution,  a  hybrid  approaeh  was 
utilized  that  combined  the  minimum  spanning  forest  and  the  modified  Clarke-Wright 
algorithms. 

LaPorte,  Nobert,  and  Taillefer  examine  three  particular  variants  of  VRP/LRPs: 
capacity  constrained  VRP,  cost  constrained  VRP,  and  cost  constrained  LRP.  They 
contend  that  producing  exact  solutions  to  these  problems  is  still  difficult,  and 
subsequently,  suggest  that  other  algorithms  are  more  applicable  in  terms  of  yielding  a 
feasible  solution.  In  particular,  they  consider  the  use  of  dynamic  programming  with 
state-space  relaxation  in  the  case  of  tightly  constrained  problems,  the  use  of  integer  linear 
programming  in  conjunction  with  constraint  relaxation  in  the  case  of  loosely  constrained 
symmetrical  problems,  and  branch  and  bound  methods  in  which  every  sub-problem  is  an 
assignment  problem  (LaPorte,  Nobert,  Taillefer,  1988).  However,  these  algorithms  can 
be  very  time  consuming  for  large  problems.  For  the  research  purposes  of  this  thesis, 
these  approaches  will  be  potentially  suitable  as  the  overall  network  structure  is  quite 
small.  Heuristic  methods,  in  contrast,  can  solve  much  larger  problems  and  often  take  into 
account  a  greater  variety  of  constraints  (LaPorte,  Nobert,  Taillefer,  1988). 

Application  to  Military  Airlift 

It  is  desired  to  examine  several  papers  that  apply  past  research  on  VRPs  based  on 
(or  that  has  potential  to  be)  a  hub-and-spoke  network  structure.  Specifically,  past 
research  conducted  on  the  optimization  of  military  airlift  assets  in  routing  cargo  is 
discussed.  Emphasis  is  placed  on  the  research’s  goals,  scenario,  variables,  and  solution 
technique. 
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Beginning  in  1994,  Naval  Postgraduate  Sehool  student  Lim  Teo-Weng  eondueted 
research  for  a  Congressional-sponsored  Mobility  Requirement  Study.  Teo-Weng 
developed  a  multi-period  Strategic  Airlift  Assets  linear  program  optimization  model 
using  the  General  Algebraic  Modeling  System  (GAMS)  (Teo-Weng,  1994).  The  overall 
goal  of  the  model  is  to  minimize  late  deliveries  and  undelivered  cargo.  As  given 
parameters,  the  demand  requirements  are  considered  to  be  known  as  are  the  number  of 
available  airlift  assets  and  airfields.  The  desired  solution  is  to  find  the  optimal 
combination  of  airlift  mission  assignments  by  number  and  type  of  aircraft  for  each  unit, 
routing  structure,  airlift  mission  start  time  and  cargo  type  to  carry  (Teo-Weng,  1994). 

The  formulated  linear  program  was  entered  in  GAMS  which  generated  a  solution  to  the 
problem. 

A  year  later,  Teo-Weng  along  with  Rosenthal  and  Morton  conducted  further 
research  into  the  strategic  airlift  assets  optimization  model  developed  by  Teo-Weng.  In 
their  paper,  they  evaluated  the  model’s  performance  over  scenarios  involving  the 
relocation  of  demand  requirements  from  Ramstein-Riyadh  to  Dhahran.  The  formulated 
model  did  not  change;  the  goal  was  to  evaluate  the  model  for  use  in  conjunction  with 
other  Air  Force  planning  tools.  The  authors  concluded  that  the  model  can  give  a 
relatively  rapid  response  to  questions  relating  to  major  mobility  issues  such  as:  1)  Are  the 
given  aircraft  and  airfield  assets  adequate  for  the  deployment  scenario?  2)  What  are  the 
impacts  of  shortfalls  in  airlift  capability?  3)  Where  are  the  system  bottlenecks  and  when 
will  they  become  noticeable?  (Morton,  Rosenthal,  Teo-Weng,  1995). 

In  1997,  RAND  briefed  the  outcome  of  a  study  requested  two  years  earlier  to 
evaluate  how  the  C-17  could  be  used  as  an  in-theater  airlifter.  Specifically,  the  study  was 
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to  determine  the  appropriate  mix  of  C- 17s  and  eivil-derived  aireraft  in  the  U.S.  military 
airlift  fleet  (Killingsworth,  Melody,  1997).  The  model  had  thousands  of  eonstraints  and 
variables.  Coneeptually,  the  formulated  model  is  fundamentally  a  VRP  altered  to  the 
speeifie  seenario  and  given  parameters.  As  in  the  researeh  eondueted  by  Teo-Weng,  the 
objeetive  of  the  model  is  to  minimize  late  deliveries  and  non-deliveries  through  a 
weighted  penalty  eoeffieient.  The  model  was  solved  using  the  GAMS  software. 

Cox  developed  a  hub-and-spoke  eombined  loeation-routing  mixed  integer 
programming  model  in  1998  as  an  alternative  to,  at  the  time,  the  U.S.  Air  Foree’s  eurrent 
direetly  delivery  methodology.  In  the  underlying  network  strueture,  he  looks  at  three 
types  of  nodes;  those  that  aet  as  hubs,  transshipment  nodes,  and  destination  nodes  (Cox, 
1998).  His  model  ineorporated  many  aspeets  from  some  of  the  authors  previously 
diseussed  in  this  literature  review,  partieularly  the  works  of  Kulkarni  and  Bhave,  and 
Baker.  The  objeetive  funetion  minimizes  lateness  and  non-deliveries,  and  utilizes  the 
following  variables: 

•  ,  is  a  binary  variable  equal  to  1  if  aireraft  with  tail  number  h  flies  the 
are  from  i  to  j 

•  Xij  is  a  binary  variable  equal  to  1  if  demand  node  j  is  supplied  by  a 
plane  based  at  depot  i 

•  Zqj  is  a  real-valued  variable,  equal  to  the  amount  of  eargo  delivered 
from  all  supply  nodes  to  transshipment  (depot)  node  j 

•  Yi  is  a  binary  variable  equal  to  1  if  any  aireraft  are  assigned  to  depot  i 

•  T* j  is  the  time  that  aireraft  h  spends  at  node  j  (Cox,  1998) 
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The  model  was  solved  using  CPLEX  3.0  whieh  generated  a  solution  in  less  than  a 
minute.  Cox  eoneluded  that  the  hub-and-spoke  model  developed  was  more  advantageous 
to  the  Air  Foree  versus  the  eurrent  method  of  eondueting  business. 

On  a  larger  seale  network,  Skipper  (2002)  analyzed  the  distribution  network  of 
U.S.  European  Command  by  formulating  a  multiple  objeetive  linear  programming  model. 
The  goal  was  to  find  an  optimal  set  of  hubs  that  would  reduee  transportation  eosts  and 
time  required  for  eargo  shipments.  Construeting  the  model  required  identifying  the 
eurrent  network  including  the  supply  node,  the  current  hub  location,  and  the  demand 
nodes.  The  variable  used  in  the  development  of  the  model  was  Xy  which  represents  the 
number  of  missions  from  node  i  to  destination  j  (Skipper,  2002).  Several  different 
scenarios  where  different  locations  acted  as  a  hub  for  the  network  were  modeled.  Excel 
was  used  as  the  solution  tool. 

That  same  year,  INFORMS  published  an  article  by  Baker,  Morton,  Rosenthal,  and 
Williams  (2002)  that  summarized  a  massive  combined  effort  on  the  part  of  the  Naval 
Postgraduate  School  and  the  RAND  Corporation  that  developed  the  NPS/RAND 
Mobility  Optimizer  (NRMO)  model.  To  date,  this  is  one  of  the  most  all-encompassing 
optimization  models  developed.  Overall,  the  model  utilizes  a  similar  set  of  general 
categories  of  variables  such  as  which  aircraft  is  delivering  to  a  particular  node  and  what 
aircraft  is  being  employed  for  making  deliveries.  However,  NRMO  breaks  these 
variables  down  into  much  more  detailed  components  resulting  in  a  significantly  higher 
number  of  variables.  GAMS  in  conjunction  with  the  parallel  CPEEX  6.0  barrier 
algorithm  was  used  to  generate  a  solution  (Baker,  Morton,  Rosenthal,  Williams,  2002). 
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Lockhart  (2005)  examined  whether  the  C-27J  aireraft  would  be  a  viable  solution 
to  help  inerease  intratheater  airlift  effieieney  in  his  researeh  paper  in  2005.  The  partieular 
methodology  utilized  in  this  researeh  represents  a  major  departure  from  the  solution 
teehniques  presented  earlier.  In  solving  the  problem,  the  researeher  eondueted  an 
analysis  of  quantitative  data  to  evaluate  potential  effieieney  improvements  if  the  C-27J  is 
added  to  U.S.  airlift  eapabilities.  Speeifieally,  a  load  faetor  ealeulation  was  developed  to 
determine  if  it  would  be  more  effieient  for  the  C-27J  to  earry  a  load  that  would  otherwise 
be  earried  on  a  minimally  loaded  C-130  (Loekhart,  2005).  Loekhart  proposes  an 
alternate  methodology  whieh  eould  be  used  to  evaluate  the  potential  replaeement  of 
several  CH-47s  with  a  eouple  of  C- 130s. 

A  more  reeent  paper  by  Rivera  (2009)  eondueted  researeh  to  determine  the  right 
mix  of  airlift  assets  that  would  be  most  suitable  for  moving  personnel  and  eargo  along  the 
Last  Taetieal  Mile.  Taking  a  more  data  analytie  approaeh,  Rivera  used  data  from  eurrent 
operations  in  Afghanistan  as  a  baseline  to  evaluate  the  effeetiveness  of  potentially 
employing  the  use  of  the  C-130J-30,  C-27J,  and  CH-47  aireraft  along  speeified  routes 
(Rivera,  2009).  His  paper  suggests  that  some  mathematieal  models  were  formulated  to 
help  form  a  solution,  but  these  are  never  presented.  However,  the  bulk  of  his  solution 
teehnique  eonsisted  of  injeeting  the  baseline  data  into  the  Air  Tasking  Effieieney  Model 
and  then  evaluating  eaeh  of  the  above  speeified  aireraft. 

Summary 

The  problem  of  redueing  wear  on  the  CH-47  blades  through  the  ineorporation  of 
additional  eargo  aireraft,  partieularly  the  C-130,  involves  many  aspeets  of  researeh  on 
hub-and-spoke  VRPs  already  aeeomplished.  The  intent  of  this  literature  review  is  to 
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demonstrate  what  mathematieal  formulations  are  eurrently  available  and  proven  to  work. 
Of  the  papers  reviewed,  it  is  evident  that  none  ean  provide  a  direet  solution  to  this 
problem.  Colleetively,  these  past  researeh  endeavors  provide  a  foundation  on  whieh  to 
formulate  a  mathematieal  integer  linear  program  unique  to  CH-47  seenario.  Although 
there  are  many  problem  formulations  for  hub-and-spoke  networks,  the  real  goal  is  to 
determine  if  some  CH-47s  ean  be  replaeed  by  C- 130s.  This  leads  to  the  VRP  as  the 
fundamental  problem  that  needs  to  be  solved  over  a  hub-and-spoke  network.  The  last 
seetion  examined  how  this  eombined  problem  has  been  applied  to  military  airlift  mobility 
problems  in  the  past  along  with  praetieal  solution  teehniques. 
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III.  Methodology 


Introduction 

It  is  desired  to  reduce  the  number  of  CH-47  hours  flown  currently  employed  in 
the  Iraq  and  Afghanistan  theater  of  operations.  In  order  to  accomplish  this,  several 
methods  were  fathomed  in  relation  to  a  mathematical  programming  approach.  An  initial 
proposed  method  was  to  begin  with  a  generalized  vehicle-routing  problem  (VRP)  and 
then  add  or  modify  constraints  specific  to  USTRANSCOM’s  scenario.  However,  VRPs 
in  the  most  general  sense  use  a  single  vehicle,  or,  if  multiple  vehicles  are  utilized,  lack 
the  detail  on  an  individual  vehicle  level  to  produce  the  information  required  for  this 
research. 

Since  this  research  focuses  on  a  multiple  aircraft-type  scenario  where  time¬ 
tracking  on  each  vehicle  is  necessary,  a  more  detailed  model  was  required  in  which 
additional  aircraft-types  could  be  added  and  easily  tracked.  This  leads  to  the  proposed 
methodology  for  this  research.  The  following  sections  detail  a  mathematical 
programming  model  with  an  underlying  hub-and-spoke  network  that  follows  a 
generalized  network  flow  problem.  By  ‘generalized,’  it  is  meant  that  an  entity  (aircraft) 
must  travel  out  of  a  source  node  (hub),  traverse  to  a  destination  or  series  of  destinations, 
and  then  proceed  to  a  sink  node  (same  hub). 

Case  Study/Network  Mapping 

The  real-world  operational  network  involves  forty-five  U.S.  bases  in  the  Iraq  and 
Afghanistan  theaters.  Due  to  the  sensitivity  of  the  data,  these  bases  will  not  be  named, 
but  a  node -base  relationship  table  can  be  established  in  any  fashion  so  long  as  those 
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nodes  acting  as  hubs  are  at  the  beginning  of  the  sequence.  These  bases  collectively  make 
up  the  entire  network,  or  node  set.  Five  of  the  bases  form  a  subset  of  nodes  that  act  as  the 
hubs  for  the  network.  The  remaining  forty  bases  are  the  destination  nodes,  or  the  nodes 
that  have  both  passenger  (PAX)  and  cargo  requirements.  Additionally,  these  destination 
nodes  can  be  formed  into  subsets  for  each  aircraft-type  as  there  are  base  restrictions  on 
whether  a  destination  node  can  accommodate  a  particular  aircraft  type.  For  this  scenario, 
the  CH-47  has  full  access  to  every  hub  and  node  in  the  network.  The  C-130,  on  the  other 
hand,  is  only  capable  of  flying  into  eighteen  nodes  from  two  hubs.  The  network  flow 
diagram  in  figure  1  captures  these  scenarios: 

Figure  12  -  General  Network  Flow 


CH-47  Hubs 
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The  top  cloud  represents  the  hubs  of  the  network  where  a  number  of  available  CH-47 
aircraft  are  stored.  Only  two  CH-47  hubs  can  accommodate  the  C-130  aircraft.  The 
bottom  cloud  consists  of  nodes  that  comprise  the  serviceable  destinations,  each  with  PAX 
and  cargo  requirements  (accounted  as  a  single  total  weight)  for  any  given  day. 

Several  assumptions  apply  to  this  network.  The  first  is  that  there  are  no  PAX  or 
cargo  requirements  between  the  five  hubs.  All  aircraft  originate  at  the  hubs  and  must 
travel  to  the  destination  nodes  to  pick  up  and  transport  passengers  and  cargo  back  to  the 
hub.  Secondly,  aircraft  do  not  travel  between  hubs.  That  is,  all  aircraft  simply  perform 
sorties  from  their  respective  hub  to  a  single  destination  or  multiple  destinations  and  back. 
Multiple  aircraft  can  service  a  single  destination.  Lastly,  there  must  be  enough  aircraft 
available  to  completely  transport  all  PAX  and  cargo. 

Integer  Program  Development 

In  accordance  with  the  prescribed  methodology,  the  variables  of  interest  are  the 
number  of  aircraft  formations  required  to  service  all  the  PAX  and  cargo  requirements  for 
a  destination  node  for  each  type  of  aircraft.  By  aircraft  formations,  it  is  meant  to  capture 
any  number  of  aircraft  that  fly  together  during  a  sortie.  For  example,  it  is  required  that 
CH-47s  always  fly  in  pairs.  Hence,  a  CH-47  formation  in  this  model  is  a  two-ship 
formation.  However,  any  number  of  aircraft  flying  together  can  be  specified  in  order  to 
model  multi-ship  formations.  The  variables  account  for  both  CH-47s  and  C-130s.  In 
order  to  establish  a  baseline  scenario  in  which  only  CH-47s  are  utilized,  the  C-130 
variables  are  ‘zeroed’  out  by  forcing  the  number  of  C-130  aircraft  available  to  zero. 
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SORTIEj  h,n,i,j 

1 1  if  arcQ,J)is  traversed  by  vehicle  type  T,  assigned  to  hub  H,  number  N 


where  < 


0,  otherwise 

fT  =  1  (CH47)  or 

2 (C130) 

H  =  1,2,.. 

,5 

N  =  1,2,.. 

,7 

I  =  1,2 . 

15 

7  =  1,2 . 

15 

additionally  I  ^  J 


NUM_REQt  =  number  of  aircraft  type  T formations  needed  to  transport  all  cargo 
requirements  from  every  destination  (*) 

Additionally,  must  be  integer  and  greater  than  or  equal  to  zero 


*Note:  NUM_REQch47  and  NUM_REQci3o  are  utilized  in  unison  so  eaeh  variable  is 
dependent  upon  the  number  of  other  vehiele  type  used. 


V_USEDt,h,n 


1  if  vehicle  type  T,  assigned  to  hub  H,  number  N  is  utilized 

0  otherwise 
H  =  1,2,..., S  if  t  =  1 

H  =  1,2  if  t  ^  2 


CE 


T,H,N,J 


amount  of  cargo  taken  from  destination  J  by  vehicle  T, 
assigned  to  hub  El,  number  N 


where 


^T,H,N,J 
must  be  >  0 
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_  (amount  of  PAX  taken  from  destination  J  by  vehicle  T, 
\  assigned  to  hub  H,  number  N 


where 


[  VT,H,N,J 

[must  be  integer  and  >  0 


A  full  scenario  will  require  291,735  variables  and  366,415  eonstraints.  This  assumes 
seven  aireraft  slots  are  available  at  eaeh  of  the  five  hubs  for  eaeh  aireraft  type,  and  at  eaeh 
of  the  45  destinations.  Again,  these  numbers  are  dependent  on  number  of  nodes,  hubs, 
aireraft  types,  and  number  of  available  aireraft. 

The  goal  for  this  mixed  integer  linear  program  (MILP)  is  to  minimize  the  total 
eosts  for  eaeh  aireraft  type  traveling  along  eaeh  hub-destination  pair.  Using  data  from 
U.S.  Transportation  Command  (USTRANSCOM)  on  flight  duration  between  hub  and 
nodes  and  DoD  estimates  on  operating  eosts  per  hour  for  eaeh  aireraft  type,  a  general  eost 
ean  be  determined  based  on  aireraft  type  and  hub-node  pair  being  traveled.  Additionally, 
sinee  the  CH-47s  fly  in  pairs,  the  eosts  assoeiated  a  CH-47  sortie  must  be  doubled.  The 
following  equation  defines  the  objeetive  funetion  for  this  ILP; 


V  V''  V''  V''  V'' 


Min  ^  Zj  Zj  Zj  *  NUM_FORMATIONt  *  TIMEtjj  *  SORTIEt,h.n.i.j) 


H  N  I 


where  I  ^  J 


The  limitations  that  oeeur  when  aireraft  type  2  (C-130)  is  utilized  are  due  to  restrietions 
on  hubs  and  destinations  that  ean  aeeommodate  the  C-130  aireraft.  In  this  seenario,  C- 
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130s  can  only  fly  out  of  hubs  one  and  two,  and  can  only  service  the  first  eighteen  nodes 
(nodes  six  through  twenty- three).  The  COSTj  vector  contains  DoD  operating  costs  per 
flight  hour  of  the  CH-47  and  C-130  in  2007.  NUM_FORMATIONt  is  a  vector  that 
specifies  the  number  of  aircraft  that  travel  together  by  type  T.  The  TYMEt;/./ parameter  is 
an  n-hy-n  matrix  where  n=  J  ox  I  and  contains  the  flight  time  (in  hours)  between  all  arc 
pairs  (T  J).  Within  the  matrix,  if  7  =  J,  then  an  arbitrarily  high  number  is  used  to  ensure 
to  no  destination  is  traveled  to  itself.  However,  these  values  are  place  holders  as 
variables  that  contain  I  =  J  pairs  are  not  created.  Cost  and  flight  time  data  is  provided  by 
USTRANSCOM. 

The  constraints  for  this  problem  have  been  broken  into  five  primary  areas.  The 
first  area  determines  the  number  of  aircraft  formations  required  to  transport  all  PAX  and 
cargo  from  all  the  destinations.  This  constraint  requires  that  there  is  sufficient  available 
aircraft  to  do  so.  For  the  PAX  and  cargo  requirements  (second  area),  only  the  total 
weight  of  both  PAX  and  cargo  is  considered  and  not  the  volume.  A  key  assumption  here 
is  that  a  “full  aircraft  capacity”  is  reached  when  fifty- five  percent  of  the  aircraft  has  been 
filled  in  terms  of  weight  (lbs).  The  third  set  restricts  aircraft  flow  throughout  the 
network.  This  set  maintains  the  network  flow  relationship  between  hub,  destination,  and 
aircraft  type.  For  example,  a  CH-47  has  access  to  every  base  that  a  C-130  does,  but  the 
C-130s  do  not  have  access  to  every  CH-47  base.  The  fourth  area  enforces  hour 
limitations  on  each  aircraft  type.  Lastly,  the  fifth  set  places  integer  restrictions  on  the 
variables,  some  of  which  are  binary. 
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Aircraft  Required. 

Three  eonstraints  determine  how  many  of  eaeh  aireraft  type  are  needed  and 
establish  whieh  of  the  available  aireraft  will  fly  sorties.  The  total  weight  eapaeities  for 
eaeh  aireraft  type  are  used  as  parameters  here  to  determine  the  number  of  aireraft 
required.  However,  sinee  this  model  does  not  aeeount  for  the  volume  spaee  of  both  eargo 
and  aireraft,  eaeh  aireraft  type’s  eapaeity  is  multiplied  by  an  adjustment  faetor  in  order  to 
aeeount  for  seenarios  where  an  aireraft’s  eargo  volume  is  fully  used,  but  is  not  at  weight 
eapaeity.  Based  on  historieal  data  from  USTRANSCOM,  this  adjustment  faetor  is  fifty- 
five  pereent. 

Constraint  la  -  Ensures  suffieient  vehiele  eapaeity  is  available; 


I 


NUM_FORMATIONt  *  WT_CAP_ADJt  *  WT_CAPt  *  NUM_REQt 


>  ^  CARGO _REQj  *  WT_PER_PAX  *  PAX_REQt 
] 


Where; 

WT_CAP_ADJt  eontains  the  adjustment  faetor  for  eaeh  vehiele  type, 
WT_CAPt  eontains  the  weight  eapaeity  by  vehiele  type, 

CARGO  REQj  is  a  veetor  speeifying  eargo  requirements  at  node  J, 

WT  PER  PAX  is  a  planning  parameter  indieating  weight  for  eaeh  passenger, 
and  PAX  REQt  is  a  veetor  speeifying  PAX  requirements  at  node  J 
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Constraint  lb  -  Ensures  total  vehicle  fleet  capacity  does  not  exceed  a  factor  for  two: 


^  NUM_FORMATIONt  *  WT_CAP_ADJt  *  WT_CAPt  *  NUM_REQt 

T 

<2*  'YjCARGO_REQj  *  WT_PER_PAX  *  PAX_REQt) 


Constraint  2a  -  Initializes  vehicles  to  be  used  based  on  need: 


V_USEDt  h,n  ^  NUM_REQt  V  T 


H  N 

Constraints  2b  -  Vehicles  used  cannot  exceed  vehicles  available  at  a  hub: 


^  V_USEDt^h,n  ^  AVAILrfl  ^T,H 

N 


where  AVAILf  ^  is  a  vector  specifying  number  of  vehicles  by  type,  assigned  to  each  hub 


PAX  and  Cargo  Constraints. 

PAX  and  cargo  are  both  accounted  for  as  a  total  weight  requirement  at  each 
destination  that  needs  to  be  transported  to  a  hub.  There  is  no  restriction  on  which  hub 
that  the  cargo  must  be  flown  to  nor  is  the  portion  of  weight  distinguished  between  PAX 
and  cargo.  The  following  constraints  restrict  an  aircraft  from  overloading  as  it  makes 
stops  at  various  destinations  as  well  as  enforces  that  the  cargo  taken  on  at  any  single 
destination  does  not  exceed  the  allotted  aircraft  weight  capacity.  This  latter  constraint  is 
likely  to  be  redundant  in  most  cases. 
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Constraints  3  -  Vehicle  cargo/PAX  weight  capacity  restriction: 

(S-T't.h.n,]  +  *  PTjhuj  j) 

J 

<  NUM_FORMATIONt  *  WT_CAP_ADJt  *  WT_CAPt 
*  V_USEDth,n  ^T,H,N 

Constraint  4  -  Vehicle  allocated  weight  capacity  for  cargo: 

^  CTt,h,n,j  ^  NUM_FORMATIONt  *  MAXjCARGO _WTr  *  V_USEDt^h,n  V  T,H,N 
J 

where  MAX_CARGO_WTt  specifies  maximum  amount  of  a  vehicle’s  weight  capacity 
allotted  to  cargo;  must  be  less  than  or  equal  to  WT_CAPt 
Constraints  5  -  Vehicle  allocated  seating  capacity  for  PAX: 

^  PTt,h,n,j  ^  NUM_FORMATIONt  *  MAX_PAXr  *  V_USEDt^h,n  V  T,H,N 
J 


where  MAX_PAXt  specifies  the  maximum  number  of  passengers  able  to  be  carried  on 

vehicle  type  T;  must  be  integer 
Constraint  6  -  Ensures  all  cargo  is  transported: 


III  CTt,h,n,j  —  CARGO _REQj  V  ] 

T  H  N 

Constraint  7  -  Ensures  all  PAX  is  transported: 


III  PTt.h.n.j  =  PAX_REQj  V  J 


T  H  N 

Constraint  8  -  Ensures  picked  up  cargo  does  not  exceed  vehicle  cargo  capacity: 

CTt,h,n,j  ^  NUM_FORMATIONt  *  WT_CAP_ADJt  *  WT_CAPt 

*  ^  SORTIEt^h,n,i,j  *  SWITCHjjj  V  T,  H,  NJ 


I 
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where  SWITCHtjj  is  a  user-defined  n-hy-n  binary  (0-1)  matrix  that  aets  as  toggle 
switehes  to  eontrol  where  a  vehiele  may  travel,  and  /  ^  J 
Constraint  9  -  Ensures  pieked  up  PAX  does  not  exeeed  vehiele  PAX  eapaeity: 

PTt,h,n,j  ^  NUM_FORMATIONt  *  MAX_PAXt 

*  ^{SORTIEt  h,n,i,j  *  SWITCHji  j  )VT,H,  N,J 

I 

where  I  ^  J 

Network  Flow  Constraints. 

Aireraft  movement  throughout  the  network  is  modeled  utilizing  a  general  fiow 
model  involving  distinet  souree  and  sink  nodes.  In  this  researeh,  a  hub  is  eonsidered  to 
be  both.  Sinee  the  network  eontains  undireeted  ares,  flows  out  of  a  node  are  eonsidered 
to  be  positive  flows  whereas  those  into  a  node  are  negative.  This  provides  a  simple  flow 
strueture  that  ensures  aireraft  leave  (if  needed)  their  assigned  hub  and  return  to  it.  In 
order  to  generate  movement,  the  hub  from  a  souree  node  perspeetive  is  given  a  flow  of  a 
single  positive  unit  and,  as  a  sink  node,  is  given  a  flow  of  a  single  negative  unit.  That  is, 
a  hub  node  has  one  unit  (an  aireraft)  that  must  be  flowed  to  itself.  Intermediate  nodes 
destinations  have  a  requirement  of  zero  flow,  whieh  forees  an  aireraft  to  stop  at  a 
destination  only  if  there  is  a  eargo  requirement.  However,  an  easily  solvable  problem 
arises  when  a  node  is  utilized  as  both  souree  and  sink,  and  the  flow  requirement  is  zero 
(instanees  when  an  aireraft  is  not  used).  If  both  souree  and  sink  network  flows  require 
zero  units  to  be  moved,  this  permits  the  intermediate  network  flows  to  potentially  transit 
an  “imaginary  aireraft”  between  an  out-flowing  and  in-flowing  are  pair,  whieh  will  take 
on  any  eargo  required  at  those  nodes  forming  the  are.  To  break  this  oeeurrenee,  the  sum 
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of  all  arc  pairs  (/,  J)  and  (/,  I)  when  I  and  J  are  greater  than  the  largest  number  of  hubs 
must  be  less  than  or  equal  to  a  unit  being  flowed,  if  any.  This  not  only  prohibits  these 
oeeurrenees,  it  also  aets  as  a  redundant  eonstraint  in  preventing  aireraft  from  returning  to 
a  previously  visited  node. 

Constraints  10a  -  Souree  node  (hub): 

^{SORTIEt^h,n,h,j  *  SWITCHt^hj)  =  V-USEDt^h,n  VT,H,N  andj 
J 

Constraints  10b  -  Destination  restrietions: 


^{-SORTIEt^h,n,i,j  *  SWITCHjjj  +  SORTIEt^h,n,j,i  *  SWITCHjj  j) 
J 

=  0  VT,H,N,J 

il 


where 


1*1 
^  H 


Constraints  lOe  -  Sink  node  (hub): 

^{-SORTIEt^h,n,j,h  *  SWITCHtj^h)  =  -y-USEDj^H.N  VT,H,N  and  J  ^  H 
J 

Constraints  lOd  -  Ensures  destinations  are  not  used  as  hubs: 

SORTIEth,n,i,j  *  SWITCHt  ij  +  SORTIEth,n,j,i  *  SWITCHjjj 
<  V_USEDt^h,n  yT.H.N.l.J 

a*i 

where  U  ^  H 

[j 
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Constraints  lOe  -  Sub-tour  elimination; 


'Y^(SORTIEt^h,n,i,j  *  SWITCHjjj  +  SORTIEt^h,n,j,i 

I 


*  SWITCHtjj) 


—  3  ,H,N ,H,K  *  ,H,N ,K,H 

*SWlTCHr^Kfl  VJ,K 


where 


rK^J 

\j^H 

U 


In  scenarios  involving  multiple  vehicle  types,  the  vehicle  type  which  utilizes  the  largest 
number  of  hubs  establishes  the  vector  length  for  the  A  VAILj  vector.  Those  vehicles  that 
do  not  utilize  a  certain  hub  are  zeroed  out  such  that  they  cannot  be  based  out  of  that 
particular  hub.  Additionally,  there  must  be  a  restriction  on  sorties  when  I  equal  J: 
Constraint  1 1  -  Ensures  no  node  self-traveling; 

SORTIEj  I  j  =  0  y  T,El,N,  I,]  and  when  I  =  J 


Aircraft  Hour  Limitations. 

Both  the  CH-47  and  the  C-130  have  total  operating  hour  restrictions.  A  six  hour 
limitation  is  placed  on  the  CH-47  while  the  C-130  can  fly  up  to  twelve  hours. 
Constraints  12  -  Vehicle  travel  duration  restrictions; 


^  ^(TIMEtj  j  X  SORTIEt^h,n,i,j)  ^  DURATION _LIMITr  V  and  I  ^  J 

I  J 

where  DURATION _LIMITt  is  a  vector  specifying  the  time  (in  hours)  limitation  for  which 

a  vehicle  type  can  operate  continuously 
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Variable  Constraints. 


The  sortie  and  aireraft  used  variables  must  be  binary  sueh  that  they  indieate 
whether  a  partieular  are  (/,  J)  was  flown  and  whether  a  partieular  aireraft  of  those 
available  was  flown.  The  variables  used  for  determining  how  many  aireraft  are  required 
are  general  integers.  Cargo  is  not  restrieted  to  be  integer,  and  inherently  assumes  that  the 
largest  pieee  of  eargo  may  be  loaded  onto  either  vehiele  types. 

Constraint  13a  -  Required  number  of  vehieles; 

NUM_REQt  >  0  and  integer  V  T 
Constraint  13b  -  Number  of  PAX  pieked  up: 

PT't,h,n,j  ^  0  integer  V  T,  H,  N,J 
Constraint  14  -  Speeifie  vehiele  utilization: 

VJJSEDt  G  [0,1]  and  integer  \/T,H,N 
Constraint  15  -  Travel  between  node  I  and  node  J: 

SORTIEt  IP  I  j  G  [0,1]  and  integer  \/ 

LINGO-based  Model  Development 

The  formulated  model  in  the  preeeding  seetion  was  solved  utilizing  LINGO  ® 
version  12.0.1.10  dated  18  January  2010.  This  formulation  eonsisted  of  16,882  integer 
variables  defining  the  number  of  aireraft  by  type  and  number  that  originate  from  a 
partieular  hub  and  travel  to  a  partieular  destination  J  from  destination  1.  This  also 
ineludes  variables  for  the  eargo  taken,  aireraft  number  refereneing,  and  the  number  of 
required  aireraft.  There  are  44,715  eonstraints  total  that  enforee  network  flow,  PAX  and 
weight  eapaeities,  and  flight  time  limitations.  The  model  was  exeeuted  on  a  HP  eomputer 
system,  model  DX5150MT  with  an  Athlon  64  X2  4400+  proeessor  and  halted  after  ten 
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minutes  of  run  time.  Initial  modeling  trials  were  eondueted  up  to  sixteen  hours  resulting 
in  no  or  minuseule  ehanges  in  the  objeetive  funetion.  In  general,  the  model  reaehes  a 
“stable”  feasible  solution  after  a  few  minutes.  Ten  minutes  of  runtime  is  seleeted  to 
ensure  an  adequate  amount  of  time  has  passed  and  no  ehanges  in  the  objeetive  funetion 
oeeurred.  On  a  side  note,  the  objeetive  funetion  ehanges  observed  while  permitting  the 
model  to  run  for  eonseeutive  hours  is  on  the  order  of  a  eouple  $1000.  In  relation  to  a 
objeetive  funetion  that  generally  totals  over  $100,000,  the  time  required  to  obtain  sueh 
small  gains  is  deemed  impraetieal. 

Through  the  objeetive  funetion  minimizing  eost  to  serviee  the  destination  bases, 
the  solution,  if  one  exists,  table  displays  the  overall  eost  for  meeting  all  destination 
requirements  and  integer  values  for  eaeh  variable.  A  striet  interpretation  of  the  sortie 
variables  should  read  as  ‘a  type  t  aireraft,  number  n  assigned  to  hub  h  flew  a  sortie  from 
destination  i  to  destination  j.’  For  praetieal  purposes,  it  may  not  be  neeessary  to  send 
multiple  aireraft  to  a  partieular  destination  even  though  a  model  solution  may  indieate 
multiple  aireraft  are  needed.  Say  for  example,  two  aireraft  are  needed  to  serviee  a 
destination  aeeording  to  this  model,  but  praetieally  only  one  aireraft  is  needed  to  perform 
two  round-trips.  In  either  ease,  this  does  not  affeet  the  overall  eost.  Ultimately,  it  is  the 
deeision  maker’s  interpretation  whether  multiple  aireraft  are  absolutely  required  or  if 
fewer  aireraft  are  required  making  multiple  trips. 

In  developing  this  model,  it  is  intended  to  eapture  all  possible  destination 
requirements  on  any  given  day.  However,  the  Iraq  and  Afghanistan  data  on  whieh  this 
model  is  based  demonstrates  that  on  an  average  day,  only  five  to  seven  destinations  have 
requirements.  This  means  that  many  of  the  eonstraints  ean  be  zeroed-out,  or  ‘switehed 


43 


off,’  greatly  reducing  the  overall  model  and  runtime  (this  can  be  a  significant  issue  when 
all  destinations  have  requirements).  Should  obtaining  a  solution  in  a  short  amount  of 
time  be  required,  the  model  does  generate  a  ‘good  enough’  (feasible)  solution  in  a  matter 
of  minutes  or  less.  Subsequently,  this  means  that  for  a  given  day,  all  constraints  will 
have  to  be  altered  for  that  day’s  destination  requirements. 

With  respect  to  updating  the  model,  some  background  on  linear  program  will  be 
required  should  it  be  desired  to  include  additional  aircraft  types  and/or  bases  that  open  or 
close  links.  However,  for  updating  the  daily  PAX  and  cargo  weight  requirements,  a 
simple  number  change  in  the  corresponding  data  table  is  sufficient.  This  also  applies  to 
updating  the  number  of  available  aircraft  stationed.  A  cautionary  note  for  updating  any 
daily  requirements  and  aircraft  availability  is  to  ensure  one  knows  the  relationship 
between  hub/destination  number  assignment  and  the  actual  base  name. 

A  last  essential  note  in  executing  this  model  is  that  it  is  mandatory  to  have  all 
daily  requirements  satisfied.  This  means  if  there  are  less  aircraft  available  than  needed  to 
meet  all  requirements,  the  model  will  not  be  able  to  generate  a  solution  (i.e.  no  feasible 
solution  will  be  found).  However,  as  noted  earlier  concerning  the  practical  use  of  aircraft 
making  roundtrips,  this  model  can  be  “short-sighted”  in  terms  of  using  few  aircraft  to 
service  all  destinations.  Still,  this  scenario,  though  a  possibility,  is  not  likely  to  occur 
given  the  current  (as  of  July  2010)  number  of  aircraft  available  at  the  nodes  and  the  low 
number  of  destinations  that  generally  have  requirements  on  a  given  day. 

Summary 

The  integer  linear  program  formulated  and  executed  in  this  research  does  have 
both  its  merits  and  limitations.  The  lack  of  connections  between  hubs  significantly 
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reduces  underlying  network  structure  and  ultimately  loses  its  similarity  with  traditional 
hub-and-spoke  problems.  This  significantly  generalizes  the  overall  model  and  permits 
relatively  simple  manipulation  through  adjusting  aircraft  availability  numbers.  Using  this 
manipulation,  an  experimental  design  can  be  developed  to  test  various  scenarios 
containing  different  aircraft  types  and  the  increase/decrease  of  aircraft  numbers  in  order 
to  evaluate  the  savings  in  CH-47  blade  hours  as  well  as  the  overall  cost.  However, 
utilizing  aircraft  with  higher  operating  costs,  such  as  the  C-130,  may  prove  to  save  CH-47 
blade  hours  at  more  monetary  expense. 

Revisiting  the  research  objectives  from  Chapter  1  on  page  2,  this  model  can 
accomplish  all  of  these.  In  regards  to  the  first  objective  of  saving  CH-47  blade  hours, 
every  aircraft  type  that  services  a  destination  results  in  airframe  savings  for  the  CH-47. 
The  model  formulated  in  this  research  generates  data  that  can  be  used  to  calculate  those 
savings.  The  second  objective  seeks  to  determine  the  amount  of  C-130  efforts  to 
generate  the  CH-47  blade  hour  savings.  Both  the  number  of  C-130s  and  the  associated 
cost  of  using  those  C-130s  yields  an  effective  determination  of  effort  required  to  obtain 
the  desired  savings.  Lastly,  the  model  itself  is  a  realization  of  achieving  the  third 
objective  of  this  research.  With  this  model,  an  optimized  mix  of  CH-47  and  C-130 
aircraft  can  be  found  while  satisfying  all  anticipated  PAX  and  cargo  requirements  in  the 
Iraq  and  Afghanistan  theater.  Additionally,  the  model  lends  itself  to  easy 
experimentation  that  permits  forced  reduction  of  CH-47  aircraft  in  order  to  generate  blade 
hour  savings. 
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IV.  Results  and  Analysis 


Introduction 

The  model  analysis  provided  in  the  following  seetion  is  based  on  results  for  a 
smaller-seale  seenario.  This  was  neeessitated  in  order  to  perform  multiple  runs  of  the 
model  and  obtain  a  feasible  solution  in  a  reasonable  amount  of  time.  Full-up  scale 
scenarios  are  conducted  as  demonstrations  of  the  model’s  ability  to  find  a  useable 
solution,  but  they  require  a  much  longer  runtime.  A  specified  scenario  where  only  CH- 
47s  are  employed  serves  as  a  baseline  for  analysis. 

Destination  node  cargo  requirements  serve  as  the  primary  drivers  in  flowing 
vehicles  throughout  the  modeled  network.  Model  results  are  based  on  achieving  100% 
cargo  movement  to  a  hub  in  the  most  cost  effective  manner  utilizing  the  air  assets 
available  at  each  hub.  An  important  note  needs  to  be  stated  in  relation  to  cargo.  From 
the  model’s  perspective,  the  number  of  PAX  in  terms  of  weight  (lbs)  and  storage  cargo 
weight  are  indistinguishable.  Cargo  requirements  at  a  destination  are  in  reference  to  the 
collective  amount  of  PAX  and  storage  cargo  weight  required  to  be  moved.  Moreover, 
where  cargo  is  moved  is  dependent  on  which  aircraft  picks  up  the  cargo. 

Scenario  Development 

The  underlying  model’s  network  is  designed  to  permit  complete  access  from  one 
node  to  any  other  node  in  the  network  aside  from  those  nodes  that  act  as  hubs.  There  is 
no  limit  prohibiting  a  user  from  entering  as  many  nodes  as  desired  so  long  as  the 
computer  hardware  contains  sufficient  memory.  This  open-ended  style  is  selected  to 
permit  this  model  to  be  employed  in  a  vast  number  of  other  scenarios  if  so  desired.  For 
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this  particular  research  scenario,  there  are  forty-five  bases  eonsidered  in  the  Iraq  theater. 
Five  of  these  bases  aet  as  major  hubs  for  the  aireraft  of  interest,  the  CH-47  and  C-130. 
With  the  forty-five  bases  alone,  2,025  ares  are  generated  that  translate  into  variables  that 
need  to  be  traeked.  However,  we  also  need  to  traek  the  different  aireraft  types,  the 
number  of  aireraft  for  eaeh  type  available,  and  to  whieh  hub  the  aireraft  is  assigned.  This 
leads  to  a  eombinatorial  variable  explosion  sueh  that  to  exhaust  every  possibility  in 
seareh  of  the  most  eost  effieient  eombination  would  be  impraetieal  in  terms  of  eomputer 
runtime. 

The  full-seale  seenario  mentioned  above  for  this  researeh  results  in  just  over 
290,000  variables.  However,  in  review  of  CH-47  flight  logs  in  the  Iraq  theater,  it  is 
found  that,  on  any  given  day,  only  six  to  seven  of  the  destination  bases  had  eargo 
requiring  movement.  To  exeeute  the  full-seale  seenario  with  those  destinations  that  had 
no  requirements  zeroed  out  would  still  be  taxing  in  terms  of  eomputer  runtime.  By 
employing  a  small-seale  seenario  with  fewer  destination  nodes,  a  feasible  solution  ean  be 
found  mueh  quieker  while  still  providing  praetieal  results  under  a  realistie  setting. 

Scenario  for  Model  Analysis. 

It  is  deeided  to  demonstrate  the  potential  savings  in  CH-47s  blade  hours  by 
utilizing  a  realistie  seenario  set-up  with  data  emulating  CH-47  and  C-130  flight  time  (in 
hours)  between  all  node  eombinations.  This  seenario  eonsists  of  fifteen  nodes  total,  the 
first  five  whieh  serve  as  major  hubs  for  the  area.  CH-47s  are  stationed  at  every  hub 
where  as  C-130s  are  only  stationed  at  the  first  two  hubs.  Additionally,  CH-47s  have 
aeeess  to  every  destination  node  in  the  network  (i.e.  nodes  six  through  fifteen).  The  C- 
130s  are  only  eapable  of  servieing  nodes  six  through  ten  due  to  airfield  restrietions  at  the 


47 


remaining  destinations.  Random  eargo  requirements  in  terms  of  thousands  of  pounds 
were  generated  for  eaeh  destination.  Aetual  aireraft  weight  eapaeities  are  implemented, 
but  redueed  by  a  55%  sealing  faetor  in  order  to  refleet  historieal  aireraft  loading  data  that 
indieates,  on  average,  that  only  55%  pereent  of  an  aireraft’s  total  weight  eapaeity  is 
utilized  due  to  an  aireraft’s  volume  limitations  being  reaehed.  Eaeh  aireraft’s  eost  per 
flying  hour  is  based  on  DoD  figures  from  2007,  data  provided  by  USTRANSCOM. 
Lastly,  the  eurrent  (as  of  2009)  number  of  CH-47s  available  at  eaeh  of  the  five  major 
hubs  in  the  Iraq  theater  are  also  implemented  here. 

The  general  seenario  established  above  is  evaluated  in  two  ways.  The  first  method 
eonsisted  of  inerementally  permitting  C-130  aireraft  to  be  available  for  use  at  partieular 
hubs,  whieh  in  effeet  allows  the  model  to  determine  if  a  C-130  should  be  utilized  if  it 
reduees  the  eost  of  eargo  delivery.  The  seeond  method  forees  a  number  of  C-130  aireraft 
to  be  utilized,  but  permits  the  model  to  determine  whieh  non-restrieted  hub  to  assign  the 
aireraft  sueh  that  eosts  ean  be  minimized.  Based  on  the  generated  eargo  requirements, 
the  aireraft  weight  eapaeities,  and  node  aeeess  restrietions,  the  most  C-130s  that  ean  be 
employed  is  three.  For  those  destinations  whieh  eannot  be  servieed  by  C-130s,  it  will 
require  a  minimum  of  eighteen  CH-47s  (or  nine  formations  of  two  CH-47s). 

Results  and  Analysis 

Scenario. 

For  the  established  seenario  settings  above,  the  two  seenarios  were  exeeuted 
several  times  eaeh  as  the  number  of  C-130  aireraft  is  ineremented.  The  first  seenario, 
whieh  referenees  the  option  of  utilizing  any  C-130  aireraft  made  available,  evaluates  eaeh 
inerement  of  C-130  at  eaeh  non-restrieted  hub  possibility.  In  the  seeond  seenario,  the 
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number  of  C-130s  forced  to  be  utilized  is  simply  evaluated  and  compared  to  the  relatable 
instance  in  scenario  one.  An  instance  where  only  CH-47s  are  flown  serves  as  the 
benchmark  with  which  to  compare  the  two  scenarios.  This  benchmark  simulates  the 
current  CH-47s  operations  in  Iraq  and  Afghanistan  theaters  such  that  no  C-130s  are  being 
flown  to  augment  the  CH-47  workload. 

For  the  benchmark  run,  no  C-130s  were  available  at  hubs  one  or  two.  The  model 
is  executed  for  ten  minutes  and  the  resulting  feasible  solution  is  evaluated.  It  is  reported 
that  fifteen  CH-47  formations  (thirty  individual  aircraft)  were  required  to  transit  all  cargo 
requirements  from  the  destinations  to  a  hub.  A  simple  calculation  of  summing  the  total 
cargo  requirement  weights  and  dividing  by  the  adjusted  weight  capacity  for  a  CH-47 
verifies  that  this  is  indeed  the  minimum  number  of  aircraft  required  to  move  all  cargo  for 
this  scenario.  Totaling  all  the  sorties  flown,  the  solution  generated  required  61.88  flying 
hours,  resulting  in  an  objective  function  cost  of  $168,499.  These  figures  are  in  reference 
to  CH-47  total  flying  hours  and  cost  for  a  day. 

With  a  captured  CH-47  workload  baseline,  the  first  scenario  is  executed  utilizing 
all  possible  C-130  increments  (up  to  three)  and  hub  assignments.  The  results  of  these 
trials  are  shown  below  in  Table  1  and  Table  2  with  reduction  values  in  relation  to  the 
baseline; 
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Table  1  -  Scenario  1  (Optional  C-130s)  Reference 

Scenario 

#of  C-130S 

#of  C-130S 

at  Hub  1 

at  Hub  2 

Sl-1 

1 

0 

Sl-2 

0 

1 

Sl-3 

2 

0 

Sl-4 

0 

2 

Sl-5 

1 

1 

Sl-6 

3 

0 

Sl-7 

0 

3 

Sl-8 

2 

1 

Sl-9 

1 

2 

Table  2-  Scenario  1  (Optional  C-130s)  Results 

Flight  Time 

Costs 

CH-47Time  Reduction 

Cost  Reduction 

1  Idl  lU 

CH-47 

C-130 

CH-47  C-130 

Total 

Time  % 

Cost 

% 

Sl-1 

49.68 

1.15 

$135,278.60 

$4,197.50 

$139,476.10 

12.2 

19.72% 

$29,023.10 

17.22% 

1.05 

Sl-2 

49.32 

1.96 

$134,298.40 

$7,154.00 

$141,452.40 

12.56 

20.30% 

$27,046.80 

16.05% 

1.05 

Sl-3 

40.62 

2.76 

$110,608.30 

$10,074.00 

$120,682.30 

21.26 

34.36% 

$47,816.90 

28.38% 

1.03 

Sl-4 

41.38 

4.95 

$112,677.70 

$18,067.50 

$130,745.20 

20.5 

33.13% 

$37,754.00 

22.41% 

1.03 

Sl-5 

43.22 

3.28 

$117,688.10 

$11,972.00 

$129,660.10 

18.66 

30.16% 

$38,839.10 

23.05% 

1.03 

Sl-6 

37.2 

3.13 

$101,295.60 

$11,424.50 

$112,720.10 

24.68 

39.88% 

$55,779.10 

33.10% 

1.14 

Sl-7 

40.02 

4.22 

$108,974.50 

$15,403.00 

$124,377.50 

21.86 

35.33% 

$44,121.70 

26.19% 

1.2 

Sl-8 

38.5 

3.08 

$104,835.50 

$11,242.00 

$116,077.50 

23.38 

37.78% 

$52,421.70 

31.11% 

1.14 

Sl-9 

40.46 

3.28 

$110,172.60 

$11,972.00 

$122,144.60 

21.42 

34.62% 

$46,354.60 

27.51% 

1.21 

Depending  on  hub  assignment,  allotting  one  C- 130  to  augment  the  CH-47’s  daily 
workload  redueed  the  workload  up  to  approximately  twenty  pereent,  or  saved  around 
twelve  hours  per  day  of  wear  on  the  CH-47’s  airframe  and  blades.  Additionally,  the 
overall  eost  for  transiting  all  eargo  requirements  is  redueed  by  sixteen  pereent.  Trials  six 
through  nine,  which  utilize  the  maximum  number  of  C-130s  that  can  be  made  available, 
demonstrate  savings  on  CH-47  blade  hours  that  average  around  thirty-six  percent. 

Overall  costs  were  reduced  by  approximately  twenty-nine  percent. 

Scenario  two,  which  required  the  designated  number  of  C- 130  aircraft  to  conduct 
sorties,  generated  similar  results  to  scenario  one  in  relation  to  the  same  increment  in  C- 
130  aircraft.  Tables  3  and  4  contain  the  results  of  scenario  two  trials; 
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Table  3-  Scenario  2  (Forced  C-130s)  Reference 

Scenario 

#of  C-130S 
required  to  fly 

S2-1 

1 

S2-2 

2 

S2-3 

3 

*S2-4 

4 

*S2-5 

5 

*S2-6 

6 

Table  4 -Scenario  2  (Forced  C- 

130s)  Results 

Scenario 

Flight  Ti 

me 

Costs 

CH-47Time  Reduction 

Cost  Reduction 

Vehideiff 

CH-47  C-130 

CH-47 

C-130 

Total 

Time  % 

Cost 

% 

S2-1 

50.54 

1.56 

$137,620.40 

$5,694.00 

$143,314.40 

11.34 

18.33% 

$25,184.80 

14.95% 

1.05 

S2-2 

43.18 

3.35 

$117,579.10 

$12,227.50 

$129,806.60 

18.7 

30.22% 

$38,692.60 

22.96% 

1.03 

S2-3 

38.26 

3.43 

$104,182.00 

$12,519.50 

$116,701.50 

23.62 

38.17% 

$51,797.70 

30.74% 

1.14 

*S2-4 

37.38 

3.67 

$101,785.70 

$13,395.50 

$115,181.20 

24.5 

39.59% 

$53,318.00 

31.64% 

1.32 

*S2-5 

37.34 

4.05 

$101,676.80 

$14,782.50 

$116,459.30 

24.54 

39.66% 

$52,039.90 

30.88% 

1.5 

*S2-6 

37.38 

4.26 

$101,785.70 

$15,549.00 

$117,334.70 

24.5 

39.59% 

$51,164.50 

30.36% 

1.67 

Similar  to  scenario  one,  the  model  is  allowed  to  mn  these  trials  for  ten  minutes 
and  the  best  feasible  solution  to  the  model  at  the  end  of  those  ten  minutes  is  reported  in 
Table  4.  Comparing  the  seenario-trial  referenees,  S2-1  in  relation  to  Sl-1  and  SI -2,  it  is 
evident  that  both  seenarios  yielded  similar  results.  Additionally,  it  should  be  observed 
that  the  model  seleeted  the  foreed  C-130  to  be  assigned  at  Hub  one  to  be  more  eost 
effieient.  Similar  observations  ean  be  made  to  S2-2  and  S2-3  in  regards  to  the  same 
number  of  C-130s  made  available  in  seenario  one.  The  asterisked  seenario  two  trials, 
four  through  six,  were  eondueted  as  additional  model  experiments.  From  a  high  level 
perspeetive,  the  use  of  subsequent  C-130  aireraft  outside  of  the  theoretieal  maximum 
number  of  these  aireraft  refleeted  relatively  no  differenee  in  the  overall  CH-47  hour 
savings  and  overall  daily  operating  eost.  Looking  at  the  results  more  in-depth,  the  exeess 
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C-130s  either  spread  the  eargo  loads  between  them  or  flew  a  sortie  along  the  shortest 
route  and  returned  to  its  hub  without  taking  on  cargo.  Figures  13  and  14  summarize  the 
CH-47  flight  time  savings  and  cost  reductions; 


Figure  13  -  CH-47  Flight  Time  Reduction 


Numberof  C-130S 


Figure  14  -  Daily  Operating  Cost  Reductions 
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Analysis  of  the  results  provided  by  the  model  demonstrates  that  signiflcant  CH-47 
blade  and  airframe  hours  can  be  saved  from  even  incorporating  a  single  C-130  aircraft. 
Even  though  the  established  scenario  here  is  small  in  comparison  to  the  total  number  of 
potential  bases  that  are  serviceable  in  the  Iraq  and  Afghanistan  theaters,  this  scenario  is 
on  par  with  the  size  of  any  given  daily  workload  network.  Applying  actual  CH-47  and  C- 
130  flight  travel  times  between  destinations  in  the  theater  network  will  still  provide 
significant  savings  in  CH-47  operating  hours. 
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Model 


Keeping  in  focus  with  the  research  objectives  of  this  research,  those  objectives  are 
revisited  in  order  to  evaluate  the  model  in  meeting  them.  Our  first  objective  is  to 
determine  how  many  CH-47  blade  hours  can  be  saved.  Based  on  scenario  and  model 
capability  demonstrations,  it  is  possible  to  determine  an  amount  of  savings  in  relation  to  a 
benchmark  involving  only  CH-47  aircraft.  The  term  “an  amount  of  savings”  is  used  to 
emphasize  that  this  model,  depending  on  platform,  has  to  be  executed  as  a  heuristic  tool 
if  the  modeled  network  is  quite  large.  From  a  practical  standpoint,  a  feasible  solution  can 
be  found  within  the  first  few  minutes  of  model  execution.  Furthermore,  a  good  estimate 
on  CH-47  blade  hour  savings  can  be  measured  with  this  model  given  that  any  aircraft  is 
forced  to  replace  a  portion  of  the  CH-47 ’s  workload.  However  these  savings  may  or  may 
not  come  at  increased  costs.  If  a  higher  capacity  yet  more  costly  aircraft  is  utilized  to 
replace  a  lighter  portion  of  the  CH-47  workload,  then  blade  hour  will  still  be  saved,  but 
the  overall  daily  costs  are  likely  to  soar.  If  a  heavier  portion  of  the  CH-47  workload  is 
replaced,  then  both  blade  hour  and  cost  savings  will  be  realized. 

The  second  research  objective  is  to  determine  what  additional  C-130  effort  is 
required  to  garner  CH-47  hour  savings.  This  objective  goes  hand-in-hand  with  the  first 
objective.  A  specific  amount  of  blade  hour  savings  is  directly  linked  to  both  the  number 
of  C-130  aircraft  utilized  and  what  sorties  those  aircraft  fly.  The  model  developed  in  this 
research  is  capable  of  generating  solutions  that  can  aid  in  determining  CH-47  savings. 
Additionally,  the  model  is  constmcted  with  sufficient  robustness  that  other  types  of 
aircraft  may  be  inserted  with  some  effort.  However,  stepping  back  and  evaluating  a 
larger  picture  of  the  CH-47  workload  problem,  what  is  not  addressed  is  the  effort 
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required  to  make  C-130s  available  from  a  joint  perspeetive.  With  that  said,  this  model 
may  be  best  utilized  as  a  high-level  planning  tool  rather  than  a  potential  daily  sortie 
planner. 

The  model  presented  and  exeeuted  during  the  eourse  of  this  researeh  is  a 
realization  of  the  third  researeh  objeetive.  At  the  start  of  this  researeh,  it  is  deeided  that  a 
method  be  eonstrueted  that  determines  the  eorreet  mix  of  C-130s  and  CH-47s  based  on 
eargo  and  passenger  requirements  and  the  number  of  C-130  eapable  and  non-C-130 
eapable  fields.  This  method  has  been  realized  in  the  form  of  an  integer-linear  program 
exeeuted  as  a  heuristie  tool.  For  large  networks,  the  best  that  this  method  ean  offer  is  a 
feasible  solution  found  within  minutes  of  model  exeeution  on  a  laptop.  Permitted 
suffieient  time  and  high-eomputing  hardware,  a  global  optimal  solution  ean  be  potentially 
found.  This  may  be  viewed  as  a  shortfall,  but  no  other  tool  exists  that  ean  yield  a 
praetieal  solution  (although  not  the  best)  and  give  the  user  signifieant  robustness  in  model 
and  network  settings. 

By  relaxing  the  variable  integer  eonstraints,  a  global  optimal  solution  ean  be 
obtained.  However,  this  solution  will  not  have  praetieal  merits  in  terms  of  variables  as 
they  will  be  fraetured.  But  we  do  gain  insight  into  a  theoretieal  lower  bound  for  the 
seenario  as  well  as  “how  good”  is  the  feasible  solution  observed  after  ten  minutes  of 
runtime.  This  theoretieal  lower  bound  must  be  equal  to  or  less  than  the  non-relaxed 
global  optimal  solution.  On  average,  eaeh  feasible  solution  obtained  from  the  seenarios 
is  within  eighty  pereent  of  the  theoretieal  optimal  solution.  Inspeetion  of  eaeh  seenario ’s 
variable  solution  report  showed  fraetured  variables,  whieh  indieates  that  the  global 
optimal  solution  found  using  relaxed  variables  is  less  than  the  theoretieal  global  optimal 
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solution  using  non-relaxed  variables.  This  suggests  that  the  obtained  feasible  solution 
may  be  even  closer  to  a  global  optimal  solution  than  eighty  percent.  Figure  15  shows  the 
obtained  feasible  solution  using  non-relaxed  variables  after  ten  minutes  of  runtime 
compared  to  a  global  optimal  solution  obtained  using  relaxed  variable  constraints. 


Figure  15  -  Bounded  Feasible  OFV 


■  Feasible  OFV  (lOmin) 

■  Relaxed  OFV  {Global  Opt) 


In  the  above  figure,  “B”  represents  the  baseline  scenario  and  “S##”  represents  the  various 
scenarios  using  the  initial  format  “S#-#.” 


55 


V.  Discussion 


Relevance  of  the  Current  Investigation 

Current  nation-building  and  stabilization  operations  in  Iraq  and  Afghanistan  have, 
to  date,  proven  to  be  the  U.S.’s  most  enduring  and  eost  intensive  efforts.  With  a 
eontinued  U.S.  presenee  and  a  sustained  need  for  supplies  and  transportation,  it  ean  be 
safe  to  antieipate  a  need  for  the  CH-47  aireraft  for  years  to  eome.  As  a  result, 
USTRANSCOM  needs  to  ensure  that  this  air  asset  is  not  exeessively  utilized  sueh  that  its 
lifespan  is  eut  short.  The  investigative  effort  eondueted  here  has  been  to  ensure  that  the 
CH-47  eontinues  to  provide  essential  airlift  eapability  as  needed  in  future  years. 

The  tool  developed  in  this  thesis  provides  essential  eapabilities  relevant  to  the 
CH-47  workloads  in  Iraq  and  Afghanistan.  Moreover,  this  tool  ean  be  applied  as  a 
heuristie  to  many  other  seenarios  involving  other  types  of  aireraft.  It  has  been 
demonstrated  to  provide  praetieal  results  whieh  ean  be  applied  to  the  various  seenarios 
involving  the  CH-47  aireraft  with  augmentation  by  the  C- 130  in  Iraq  and  Afghanistan. 
Reflections  on  Obtained  USTRANSCOM  Data 

Data  provided  by  USTRANSCOM  provided  signifieant  insight  into  the  CH-47 
operations  in  Iraq  for  2009  and  2010.  Mueh  of  the  underlying  network  is  derived  based 
on  where  CH-47s  eondueted  sorties.  In  order  to  avoid  olassifieation,  notimal  nodes  with 
emulated  aireraft  travel  times  are  utilized.  To  replieate  the  underlying  CH-47  networks  in 
Iraq,  one  would  simply  need  to  assign  real  air  base  names  to  nodes  and  enter  aetual 
aireraft  flight  times  in  the  model’s  data  matriees. 
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Since  it  was  decided  to  construct  a  general  model  capable  of  capturing  the  CH-47 
networks  in  Iraq  and  Afghanistan,  no  data  is  actually  used  to  generate  the  network(s). 

The  data  required  to  construct  exact  CH-47  networks  is  inputted  by  the  user,  effectively 
mitigating  the  need  to  evaluate  data  integrity  from  the  perspective  of  this  research. 
However,  this  could  become  an  issue  for  the  user.  The  model  generates  the  required 
variables  and  constraints  based  on  the  data  entered  by  the  user. 

In  essence,  the  data  from  USTRANSCOM  is  invaluable  in  determining  the 
model’s  aircraft  engineering  parameters,  flying  hour  costs,  and  network  flow  and 
accessibility.  The  data  contained  records  of  CH-47  flight  sorties  along  with  cargo  and 
PAX  taken  on  or  dropped  off.  These  records  provided  a  great  deal  of  insight  into  the 
utilization  of  the  aircraft  in  the  Iraq  theater.  This  utilizations  played  a  factor  in 
determining  how  best  to  adjust  the  weight  capacity  of  the  aircraft  to  account  for  volume 
restrictions. 

Perspective 

From  an  employment  perspective,  the  model  is  intended  to  minimize  the  cost  in 
an  established  network  given  available  vehicle  resources  while  transiting  ah  node 
requirements.  At  the  most  basic  level,  the  mode  attempts  to  find  the  best  mix  of  vehicles 
given  capacity  and  time  constraints,  and  where  those  vehicles  should  travel.  It  is  possible 
to  utilize  this  model  as  a  daily  planning  tool  as  it  yields  a  feasible  solution  rapidly. 

Evaluating  the  developed  model  from  a  mathematical  perspective,  it  combines 
aspects  of  minimum  cost  flow  and  vehicle  routing  models  on  a  hub-and-spoke  network 
uniquely  tailored  to  provide  balanced  solutions  for  the  CH-47  workload.  However, 
despite  being  designed  with  variables  and  parameters  for  the  CH-47  and  C-130  aircraft. 
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the  model  ean  be  easily  manipulated  to  represent  any  kind  of  vehiele  flow.  A  major 
drawbaek  is  the  number  of  variables  required  to  traek  vehiele  designations  and  route 
eombinations  that  ean  rapidly  explode. 

On  a  researeher  note,  model  development  took  several  iterations  and  a  signifieant 
amount  of  time  to  develop.  The  lengthy  endeavor  is  primarily  due  to  attempts  to  keep  the 
model  open-ended  enough  sueh  that  it  is  not  limited  to  the  exaet  CH-47  and  C-130 
networks  eurrently  employed  in  Iraq.  Similarly,  it  is  not  desired  to  limit  the  model  to 
those  partieular  aireraft.  The  end  result  of  this  researeh  is  a  general  model  in  whieh  the 
user  is  free  to  set  their  parameters,  vehieles,  and  network.  The  only  restrietions  to  the 
model  in  regards  to  being  open-ended  are  that  eross-flow  between  hubs  is  non-existent 
and  all  eargo  requirements  at  the  destinations  must  be  transported  to  a  hub. 

Conclusion 

Replaeing  any  portion  of  the  CH-47  workload  with  a  larger  eapaeity  aireraft  will, 
intuitively,  net  savings  on  the  CH-47  airframe  and  blade  hours.  The  developed  model 
links  the  CH-47  workload  replaeement  through  the  use  of  other  aireraft  types  and 
generates  a  praetieal  solution.  Through  investigative  studies  of  network  flow  models, 
vehiele  routing,  and  hub-and-spoke  networks,  this  researeh  has  tied  the  areas  together 
into  a  single  mathematieal  program  designed  to  reduee  the  CH-47  workload.  This  study 
eoneludes  by  offering  an  integer-linear  program  that  ean  be  used  as  a  heuristie  that 
attempts  to  diseover  the  minimal  amount  of  aireraft  use  required  to  move  100%  of  a  set 
of  eargo  requirements,  given  a  set  of  hubs  and  destinations  that  inelude  network 
restrietion. 
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Appendix  A,  LINGO-Based  Model 


MODEL: 

TITLE  Optimal  CH-47  and  C-130  Workload  Balance; 

I  -k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

! ***Variables***; 

!  SORTIE (T, H, N, I, J)  =  1  if  a  vehicle  formation  of  type  T,  number  N 

stationed  at  hub  H  travels  from  node  I  to 
node  J; 

!  0  otherwise; 

!  where  T  =  1,  2,  'user  specified'; 

!  N  =  1,  2,  'user  specified'; 

!  H  =  1,  2,  'user  specified'  and 

less  than  [max(I,J)  -  1]; 

!  I,  J  =  1,  2,  'user  specified'; 

!  CT{T,H,N,J)  =  amount  of  cargo  taken  from  node  J  by  a  vehicle  of 

type  T,  number  N,  stationed  at  hub  H; 

!NOTE:  CT  refers  to  'Cargo  Taken'  and  cargo  units  are  user-defined  in 
the  vehicle  settings; 

!  PT{T,H,N,J)  =  number  of  PAX  taken  from  node  J  by  a  vehicle  of 

type  T,  number  N,  stationed  at  hub  H; 

!NOTE:  PT  refers  to  'PAX  Taken'; 

!  NUM_REQ{T)  =  the  minimum  number  of  type  T  vehicle  formations 

required  to  move  all  PAX/cargo  requirements; 

!  V_USED (T, H, N)  =  1  if  a  vehicle  formation  of  type  T,  number  N  was 

used  to  conduct  sorties  out  of  hub  H; 

!  0  otherwise; 

I********************************************************************* 

SETS: 

DESTINATION:  CARGO_REQ,  PAX_REQ; 

HUB:  AVAIL_CH47,  AVAIL_CI30; 

VEHICLE_TYPE :  NUM_HUBS,  COST,  WT_CAP,  MAX_CARGO_WT ,  MAX_PAX, 

NUM_FORMATION,  WT_CAP_ADJ,  DURATION_LIMIT,  NUM_REQ 
TOTAL_VEHICLE_TIME,  TOTAL_VEHICLE_COST , 

T  0  T AL_C ARGO_MO VE  D ,  T 0 T AL_P AX_MO VE  D ; 

VEHICLE_NUM; 

DxD (DESTINATION,  DESTINATION):  TIME_CH47,  TIME_CI30,  SWITCH_CH47 

SWITCH_CI30; 

TxH {VEHICLE_TYPE,  HUB):  AVAIL; 

HxN{HUB,  VEHICLE_NUM) ; 

TxHxN {VEHICLE_TYPE,  HUB,  VEHICLE_NUM) :  V_USED; 

TxDxD {VEHICLE_TYPE,  DESTINATION,  DESTINATION) :  TIME,  SWITCH; 
HxNxD{HUB,  VEHICLE_NUM,  DESTINATION); 

TxHxNxD {VEHICLE_TYPE,  HUB,  VEHICLE_NUM,  DESTINATION):  CT,  PT; 
HxNxDxD{HUB,  VEHICLE  NUM,  DESTINATION,  DESTINATION); 


59 


TxHxNxDxD {VEHICLE_TYPE,  HUB,  VEHICLE_NUM,  DESTINATION, 
DESTINATION) :  SORTIE; 

ENDSETS 

I  -k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k* 

•  r 

DATA: 

!***HUB  AND  DESTINATION  SETTINGS***; 

DESTINATION  =  NI  N2  N3  N4  N5  N6  N7  N8  N9  NIO  Nil  NI2  NI3  NI4  NI5; 

ISets  total  number  of  nodes  (combined  hubs  and  destinations)  for  the 
network; 

HUB  =  Hi  H2  H3  H4  H5; 

ISets  the  maximum  number  of  hubs  for  the  network; 

!***CARGO  AND  PAX  SETTINGS***; 

CARGO_REQ  =  @FILE {' Cargo . PRN ') ; 

IThis  file  contains  all  the  cargo  requirements  (in  1000  lbs)  for  each 
destination; 

PAX_REQ  =  @FILE { ' PAX. PRN' ) ; 

IThis  file  contains  all  the  PAX  requirements  (in  1000  lbs)  for  each 
destination; 

!***VEHICLE  DATA***; 

I ******gH47 ******. 

TIME_CH47  =  @FILE ( 'CH47_Times.PRN' ) ; 

! These  are  the  times  in  hours  req'ed  for  a  CH-47  to  fly  from  node  I  to 
node  J; 

SWITCH_CH47  =  @FILE ( ' CH47_Switch_Matrix . PRN ' ) ; 

IThis  binary  matrix  controls  where  a  CH-47  can  fly; 

AVAIL_CH47  =  @ FILE (' CH4 7_Availabi lity . PRN ') ; 

IThis  file  contains  the  number  of  CH-47  stationed  at  each  node; 

I ******C430******; 

TIME_C130  =  @FILE ( 'C130_Times.PRN' ) ; 

! These  are  the  times  in  hours  req'ed  for  a  C-130  to  fly  from  node  I  to 
node  J; 

SWITCH_C130  =  @FILE ( 'C130_Switch_Matrix.PRN' ) ; 

IThis  binary  matrix  controls  where  a  C-130  can  fly; 

AVAIL_C130  =  @FILE ( 'C130_Availability.PRN' ) ; 

IThis  file  contains  the  number  of  C-130  stationed  at  each  node; 

!***VEHICLE  SETTINGS***; 

VEHICLE_TYPE  =  CH47  CI30; 

! Enter  a  vehicle  type  ID  designator  here; 

VEHICLE_NUM  =  A1  A2  A3  A4  A5  A6  A7 ; 

INOTE:  This  set  must  contain  a  number  of  vehicle  slots  equal  or  greater 
than  the  highest  number  of  required  vehicle  for  either  type, 
i.e.  {Al, A2, . . . , AN}  where  N  =  MAX (NUM_REQ (Tl) ,  NUM_REQ (T2 ),...) ; 
NUM_HUBS  =52; 

! Establishes  the  number  of  hubs  in  the  network  starting  at  node  1,  node 
2,  . . .,  node  J  for  each  vehicle  type; 

COST  =  2723  3650; 

! Established  the  cost  of  vehicle  operation  per  hour  by  type; 

WT_CAP  =  8.5  45.0; 

ISets  the  weight  (in  1000  lbs)  capacity  for  the  CH-47  VEHICLE; 

MAX  CARGO  WT  =  8.5  45.0; 
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!Sets  a  maximum  allocated  amount  of  weight  for  cargo  by  vehicle  type, 
must  be  equal  or  less  than  WT_CAP; 

MAX_PAX  =  8  150; 

!This  is  the  maximum  number  of  PAX  that  can  be  carried  by  vehicle  type; 
NUM_FORMATION  =21; 

! Establishes  the  number  of  vehicles  that  travel  together  in  a  formation 
or  group  by  type; 

WT_CAP_ADJ  =  0.55  0.55; 

! Adjusts  a  vehicle's  weight  capacity  to  account  for  the  vehicle's 
volume  limits  by  type; 

DURATION_LIMIT  =  6  12; 

!Sets  the  duration  limit  on  how  long  a  vehicle  can  continuously 
operate; 

WT_PER_PAX  =  0.4; 

!Sets  the  allocated  weight  (in  1000  lbs)  for  a  single  PAX; 

ENDDATA 

!NUM_REQ{1)  =  17; 

! NUM_REQ { 2 )  =  0  ; 

! These  settings  force  a  number  of  specific  vehicle  type  to  be  used, 
must  be  equal  or  less  than  the  total  number  of  a  vehicle  type 
available; 


•  r 

! ***PREPARATION  OPERATIONS***; 

lUsed  to  read  in  vehicle  travel  time  data; 


@FOR{TxDxD{T,I,  J) 

IT 

#EQ# 

1: 

TIME  (T,  I,  J) 

= 

TIME 

_CH47 (I, J) ) ; 

@FOR{TxDxD{T,I,  J) 

IT 

#EQ# 

2  : 

TIME (T, I, J) 

= 

TIME 

_C13  0  (I,  J)  )  ; 

lUsed  to  read  in 

vehicle 

node  access 

@FOR{TxDxD{T, I, J) 

IT 

#EQ# 

1: 

SWITCH (T, I, 

J) 

=  SWITCH  CH47 (I, 

@FOR{TxDxD{T, I, J) 

IT 

#EQ# 

2  : 

SWITCH (T, I, 

J) 

=  SWITCH  C130 (I, 

lUsed  to  read  in  vehicle  availability  data; 

@FOR{TxH{T,H) IT  #EQ#  1: 

AVAIL{T,H)  =  AVAIL_CH47 (H) ) ; 

@FOR{TxH{T,H) IT  #EQ#  2: 

AVAIL{T,H)  =  AVAIL_C130 (H) ) ; 

I*********************************************************************, 

•  t 

! ***OBJECTIVE  FUNCTION***; 

MIN  =  @SUM{TxHxNxDxD{T,H,N, I, J)  I  I  #NE#  J: 

COST{T)  *  NUM  FORMATION{T)  *  TIME{T,I,J)  *  SORTIE { T , H, N, I ,J)); 


I*********************************************************************, 

•  t 

!  ***CONSTRAINTS***; 

!  ******Vehicle  availability  and  requirement  constraints’ 

!  *********Constraint  la*********,- 


.******. 
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@SUM{VEHICLE_TYPE (T) : 

NUM_FORMATION (T)  *  WT_CAP_AD J { T )  *  WT_CAP{T)  *  NUM_REQ (T) )  > 

@SUM (DESTINATION (J) : 

CARGO_REQ{J)  +  WT_PER_PAX  *  PAX_REQ { J) ) ; 

I *********Constraint  Ib *********; 

@SUM{VEHICLE_TYPE (T) : 

NUM_FORMATION (T)  *  WT_CAP_AD J { T )  *  WT_CAP{T)  *  NUM_REQ (T) )  < 

2  *  @ SUM (DESTINATION (J) ; 

CARGO_REQ(J)  +  WT_PER_PAX  *  PAX_REQ(J)); 

I  ******** *Constraint  2a*********; 

@FOR(VEHICLE_TYPE (T) ; 

@SUM(HxN (H,N) :  V_USED (T, H, N) )  =  NUM_REQ(T)); 

I *********Constraint  2b *********; 

@FOR(TxH (T,H) : 

@SUM(VEHICLE_NUM(N) :  V_USED ( T , H , N) )  <=AVAIL(T,H)); 

! ******PAX/cargo  vehicle  capacity  constraints******; 

I  *********Constraint  3*  ******** ; 

@FOR(TxHxN(T,H,N) : 

@SUM (DESTINATION (J) : 

CT(T,H,N,J)  +  WT_PER_PAX  *  PT(T,H,N,J))  <= 
NUM_FORMATION (T)  *  WT_CAP_AD J ( T )  *  WT_CAP(T)  * 

V_USED (T, H, N) ) ; 

I  ******* **Constraint  4*********; 

@FOR(TxHxN(T,H,N) : 

@SUM (DESTINATION (J) : 

CT (T, H, N, J) )  <= 

NUM_FORMATION (T)  *  MAX_CARGO_WT (T)  *  V_USED ( T , H, N) ) ; 

I *********Constraint  5*********; 

@FOR(TxHxN(T,H,N) : 

@SUM (DESTINATION (J) : 

PT  (T,  H,  N,  J)  )  <= 

NUM_FORMATION (T)  *  MAX_PAX(T)  *  V_USED ( T , H, N) ) ; 

I *********Constraint  6*********; 

@FOR (DESTINATION (J) ; 

@SUM(TxHxN(T,H,N) :  CT(T,H,N,J))  =  CARGO_REQ ( J) ) ; 

I *********Constraint  7*********; 

@FOR (DESTINATION (J) : 

@SUM(TxHxN(T,H,N) :  PT(T,H,N,J))  =  PAX_REQ(J)); 

! *********Constraint  8*********; 

@FOR(TxHxNxD(T,H,N, J) : 

CT(T,H,N,J)  <= 

NUM_FORMATION (T)  *  WT_CAP_AD J ( T )  *  WT_CAP(T)  * 

@SUM (DESTINATION (I)  I  I  #NE#  J: 

SORTIE (T,H, N, I, J) ) ) ; 
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1  *********Constraint  9*********; 

@FOR{TxHxNxD{T,H,N,  J)  : 

PT{T,H,N,J)  <= 

NUM_FORMATION (T)  *  MAX_PAX{T)  * 

@SUM (DESTINATION (I)  I  I  #NE#  J: 

SORTIE (T, H, N, I, J) ) ) ; 

I  ******Network;  constraints******; 

!  *********Constraints  lOa,  lOb,  lOc,  lOd,  lOe,  and  lOf*********; 
@FOR{TxHxN(T,H,N)  : 

@SUM (DESTINATION (J) |J  #NE#  H: 

SORTIE(T,H,N,H, J) )  =  V_USED(T,H,N); 

@FOR (DESTINATION (J) |J  #NE#  H: 

@SUM(DESTINATION (I)  I  I  #NE#  J;  -SORTIE (T, H, N, I , J)  )  + 
@SUM(DESTINATION (I)  I  I  #NE#  J;  SORT IE ( T , H, N, J, I )  )  =  0); 

@SUM (DESTINATION (J) |J  #NE#  H: 

-SORTIE (T,H,N,J,H) )  =  -V_USED(T,H,N); 

@FOR(DxD (I, J)  I  I  #NE#  J  #AND#  I  #NE#  H  #AND#  J  #NE#  H: 

SORTIE (T,H, N, I, J)  +  SORTIE ( T , H, N, J, I )  <=  V_USED ( T , H, N) ) 

@FOR(DxD (K, J) |K  #NE#  H  #AND#  J  #NE#  H: 

@SUM(DESTINATION(I) |I  #NE#  H  #AND#  I  #NE#  J; 

SORTIE (T, H,N, I, J)  +  SORTIE ( T , H, N, J, I ) )  <= 

3  -  SORTIE (T,H,N,H,K)  -  SORTIE(T,H,N,K,H)); 

@FOR(DxD (I, J)  I  I  #NE#  J; 

SORTIE (T, H,N, I, J)  <=  SWITCH ( T , I , J)  ))  ; 

I  ******* **Constraint  ll*********,- 
@FOR(TxHxNxDxD (T, H, N, I, J)  I  I  #EQ#  J: 

SORTIE (T, H,N, I, J)  =  0) ; 


!******Hour  constraints******; 

!  ******* **Constraint  12*********,- 
@FOR(TxHxN(T,H,N) : 

@SUM(DxD (I, J) |J  #NE#  I: 

TIME(T,I,J)  *  SORTIE (T, H, N, I, J)  *  SWITCH ( T , I , J) )  <= 

DURATION_LIMIT (T)  )  ; 

I ******variable  constraints******; 

1  *********Constraint  13*********,- 
@FOR(VEHICLE_TYPE (T) : 

@GIN (NUM_REQ (T) ) ) ; 

@FOR(TxHxNxD(T,H,N, J) : 

@GIN(PT(T,H,N, J) ) ) ; 

1  *********Constraint  14*********,- 
@FOR(TxHxN(T,H,N) : 

@BIN(V  USED(T,H,N)  )  )  ; 
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1  *********Constraint  15*********,- 
@FOR{TxHxNxDxD{T,H,N,I, J) |J  #NE#  I: 

@BIN (SORTIE {T,H,N, I, J) ) ) ; 

I  ic-kicic-kic-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

!***ADD'L  OUTPUT  CALCULATIONS***; 

INOTE:  To  reduce  model  runtime,  disable  (comment-out)  the  below 
calculations ; 

! Total  weight-carrying  capacity  of  entire  vehicle  fleet; 

AVAIL_WT_CAP  =  @SUM(TxH (T,H) : 

NUM_FORMATION (T)  *  WT_CAP_AD J ( T )  *  WT_CAP(T)  * 

AVAIL (T, H) ) ; 

! Weight-carrying  capacity  used  for  the  established  scenario  (what  the 
model  actually  used) ; 

USED_WT_CAP  =  @SUM(VEHICLE_TYPE (T) : 

NUM_FORMATION (T)  *  WT_CAP_AD J ( T )  *  WT_CAP(T)  * 

NUM_REQ (T) ) ; 

! Weight-carrying  capacity  required  for  the  established  scenario  (what 
is  actually  needed) ; 

REQ_WT_CAP  =  @SUM (DESTINATION (J) : 

CARGO_REQ(J)  +  WT_PER_PAX  *  PAX_REQ(J)); 

IVehicle  use  efficiency  (gauges  model  efficiency  -  the  closer  this 
ratio  is  to  1.0,  the  more  efficient); 

VEHICLE_EFF  =  USED_WT_CAP  /  REQ_WT_CAP; 

! Total  vehicle  travel  time  based  on  sorties  conducted; 
@FOR(VEHICLE_TYPE (T) ; 

TOTAL_VEHICLE_TIME (T)  = 

@SUM(HxNxDxD (H,N, I, J) I J  #NE#  I: 

NUM_FORMATION (T)  *  TIME(T,I,J)  *  SORTIE (T, H, N, I , J) )) ; 

! Total  vehicle  costs  based  on  sorties  conducted; 

@FOR(VEHICLE_TYPE (T) : 

TOTAL_VEHICLE_COST (T)  = 

@SUM(HxNxDxD (H,N, I, J) I J  #NE#  I: 

COST(T)  *  NUM_FORMATION (T)  *  TIME(T,I,J)  * 

SORTIE (T,H, N, I,  J)  )  )  ; 

! Total  cargo  moved  by  vehicle  type; 

@FOR(VEHICLE_TYPE (T)  : 

TOTAL_CARGO_MOVED (T)  =  @SUM (HxNxD (H, N, J) :  CT ( T , H, N, J) ) ) ; 

! Total  number  of  PAX  moved  by  vehicle  type; 

@FOR(VEHICLE_TYPE (T) : 

TOTAL_PAX_MOVED (T)  =  @ SUM (HxNxD ( H, N, J) :  PT ( T , H , N, J) ) ) ; 

I********************************************************************* 

END 
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Appendix  B,  Blue  Dart 


How  can  we  reduce  the  CH-47  workload  in  the  Iraq  AOR  while  maintaining 
fiscal  responsibility?  The  answer  is  to  replace  the  heaviest  portions  of  the  CH-47 
workload  with  few  high-capacity  airlift  assets.  Including  even  a  single  high-capacity 
aircraft,  despite  incurring  steeper  operation  costs,  is  significantly  more  efficient  than  the 
collective  operation  costs  of  the  CH-47s  replaced. 

Current  Iraq  and  Afghanistan  internal  supply  and  delivery  operations  are  handled 
largely  by  the  CH-47  aircraft.  These  operations  consist  of  movement  of  cargo  and 
passengers  out  of  centralized  hubs  to  various  forward  operating  bases,  and  vice-versa.  A 
subset  of  hubs  and  FOBs  include  overlapping  operations  performed  by  other  airlift  assets, 
namely  the  C-130.  Through  optimization  modeling,  significant  overlap  can  be  reduced 
by  transferring  a  portion  of  the  CH-47  workload  to  the  C-130.  This  ultimately  will 
increase  the  CH-47  lifespan  and  reduce  daily  operating  costs. 

Intuitively,  replacing  any  portion  of  the  CH-47 ’s  daily  workload  will  net  savings 
on  the  airframe,  blade  hours,  and  maintenance.  However,  this  is  will  not  necessarily 
translate  into  reduction  of  daily  operation  costs.  We  can  better  ensure  costs  are  trimmed 
by  targeting  the  heavy  workload  portions.  An  optimization  model  helps  us  find  these 
targeted  workloads. 

A  mixed-integer  linear  program  has  been  developed  to  find  a  balance  in  the  CH- 
47  workload  with  C-130  augmentation  while  minimizing  daily  operating  costs.  This 
model  permits  ‘user-defined’  network  structures  with  which  to  evaluate  costs  utilizing  a 
set  of  airlift  assets.  Currently,  the  model  captures  “pick-up  and  return  to  hub”  cargo  and 
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passenger  seenarios.  This  eapability  mimies  eighty  pereent  of  eurrent  CH-47  operations 
in  Iraq  and  Afghanistan. 

The  CH-47  seenario  in  Iraq  has  been  implemented  in  the  model  to  demonstrate 
the  potential  value  to  the  U.S.  military.  In  this  seenario,  only  a  few  hubs  and  FOBs  have 
the  eapability  to  aeeommodate  the  C-130  aireraft.  This  implieation  suggests  that  the  most 
CH-47  workload  replaeed  by  the  C-130  eonsist  of  loeations  that  are  “C-130  eapable.” 
Several  seenarios  were  modeled,  with  eaeh  seenario  relating  to  a  higher  ineremented 
number  of  C-130s  utilized.  The  savings  are  rather  signifieant.  Based  on  simulated  eargo 
and  passenger  requirements  for  a  given  day,  potentially  between  twenty-six  and  forty- 
three  pereent  in  CH-47  flight  time  ean  be  saved  in  relation  to  the  time  that  would  be 
required  if  operations  were  strietly  handled  by  the  CH-47.  Similarly,  daily  operating 
eosts  ean  potentially  be  redueed  between  twenty-five  and  thirty-six  pereent. 

The  results  of  this  simulated  seenario  in  Iraq  attest  to  the  potential  value  that 
optimization  modeling  ean  provide  when  applied  to  operations  eontaining  a  high  number 
of  low  eapaeity  airlift  assets.  As  the  eurrent  U.S.  nation-building  and  stabilization 
endeavors  in  Iraq  and  Afghanistan  are  sure  to  be  long,  the  day-to-day  savings  garnered  by 
the  MILP  employed  will  prove  signifieant  for  the  duration.  On  a  similar  note,  the  U.S. 
will  need  to  ensure  that  the  lifespan  of  its  eurrent  airlift  assets  are  lengthened. 
Additionally,  suppose  one  desires  to  manage  a  large  number  of  high  eapaeity  airlift  assets 
more  effieiently.  Optimization  modeling  ean  yield  similar  results  as  well  by  replaeing  the 
lighter  workload  portions  often  eondueted  by  high  eapaeity  aireraft  with  a  more  tailored 
eapaeity  airlift  solution. 
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The  bottom  line  is;  we  must  eontinue  to  meet  our  mission  requirements  under 
eontinually  tightening  fiseal  eonstraints.  Targeting  the  appropriate  workload  densities 
with  the  eorreetly  tailored  airlift  assets  will  ensure  eontinued  mission  sueeess  and  out 
down  on  daily  oosts.  Optimization  modeling  ensures  we  have  the  right  mix  of  airlift 
assets  for  the  mission,  operate  at  lower  oosts,  and  preserves  are  airframes. 


The  views  expressed  in  this  article  are  those  of  the  author  and  do  not  reflect  the  official 
policy  or  position  of  the  United  States  Air  Force,  Department  of  Defense,  or  the  US 
Government. 
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Appendix  C.  Storyboard 


Optimal  CH-47  and  C-130  Workload 
Balance 


INTRODUCTION 

Since  (he  U.S.  began  its  nation-building  and 
stabilization  endeavors  in  Iraq  and  Afghanistan,  it 
has  heavily  relied  on  the  CH^7  to  move  cargo  and 
passengers  to  and  from  centralized  hubs  and 
forward  operating  bases.  The  austere  terrain  and 
high  OPS  tentpo  at  these  locations  have  placed  a 
heavy  burden  and  rapid  wear  on  the  CH^7  airframe. 
In  order  to  expand  the  lifespan  and  reduce 
maintenance  operations  of  this  airframe. 
USTRANSCOM  required  a  study  to  evaluate  the 
potential  CH47  airframe  and  cost  savings  by 
augmenting  its  workload  with  the  C-130.  The  set  of 
hubs  and  destinations  which  (he  CH-47  service  are 
not  all  accessible  by  the  C-130.  As  such,  the  portion 
of  workload  to  be  augmented  encompasses  those 
hubs  and  destinations  which  are  capable  of 
receivingC-130$.  Eighty  percent  of  the  cargo  and 
passenger  requirements  are  between  the  various 
hubs  artd  rlestiitations.  The  remaining  twenty 
percent  involves  movement  of  requirements 
between  the  destinations.  It  is  desired  to  evaluate 
the  trade-off  between  flight  time  hours  saved  on  the 
CH-47  airframe  and  the  costs  for  incorporating  the 
C-130.  Replacing  high-density  CH-47  workload 
portions  that  require  multiple  aircraR  with  few.  yet 
more  costly  C-130s  will  ultimately  reduce  CH-47 
wearand  reduce  total  cost. 


RESEARCH  GOAL 

The  goal  of  this  research  is  to 
reduce  wear  on  the  CH-47 
blades  and  airframe  by  finding 
an  optimal  workload  balance 
between  the  CH-47  and  the 
inclusion  of  the  C-130.  This 
research  determines  the  mix  of 
optimal  airlift  assets  for  Iraq. 


Captain  Dustin  P.  Jones 
Advisor:  Dr.  James  T.  Moore 
Department  of  Operational  Sciences  (ENS) 
Air  Force  Institute  of  Technology 

MODEL  FRAMEWORK 


MIXED-INTEGER  LINEAR  PROGRAM 


It  network  flow. 
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Objactiva  Function: 


The  objective  function  yields  a  total 
daily  operating  costs  based  on  the 
set  of  vehicles  utilized: 

-  Cost  per  operating  hour 


-  Time  to  traverse  arc(/,J) 
Constraints  include  five  areas 
aircraft  availability.  cargo/PAX 


Scenario  Evaluated 

The  CH-47  network  in  Iraq  is  (he  scenario 
implemented  in  this  research  used  to 
evaluate  the  mixed-integer  linear  program 
model  developed.  Based  on  the  number  of 
C-130  accessible  hubs  and  destinations,  at 
most  three  C-130s  can  be  incorporated  to 
reduce  CH-47  flight  time. 

Full  Combinatorial  Testing 


Inclusion  of  the  C-130s  demonstrate  significant 
flight  time  savings  on  the  CH-47  as  well  as 
reduction  in  daily  operaling  costs. 


RESULTS  AND  CONCLUSIONS 

Augmenting  current  CH-47  operations  in  Iraq 

and  Afghanistan  with  C-130s  will  have 

substantial  payoff: 

•  An  optimal  blend  of  CH-47  and  C-13D  airlift 
assets  reduces  up  to  37%  flight  time  wear  on 
the  CH-47  on  a  day-to-day  basis 

•  TargetingtheheavyCH-47workloadportions 
with  C-130  replacement  can  reduce  daily 
operating  costs  by  nearly  30% 

Future  Research; 

•  Augmentmodel  with  additionatairframes 

•  Capture  inter-destination  nodes  cargo/PAX 
requirements 

•  Include  time-windows  to  account  for  longer 
missions  and  operational  restrictions 


Sponsor 

U.S  Transpon  jmmand  (USTRANSCOM) 
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